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ABSTRACT 

Three  different  approaches  have  been  developed  by  groups  headed 
by  Hansch,  Holmes  and  by. Free  and  Wilson  to  the  problem  of  quantitative 
structure-activity  relationships  (QSAR)  of  drugs  to  biological  processes 
and  these  have  been  reviewed  and  the  strengths  and  weaknesses  cf  each 
method  annotated.  — 

The  five  general  equations  developed  by  Hansch  from  extra- 
thermodynamic  considerations  are  presented  as  well  as  the  ones  supple- 
mented by  inclusion  of  steric  substituent  constants  and  dunu\y  or  indicator 
parameters.  Examples  of  10  representative  type  equations  are  presente, ' 
for  the  calculation  of  biological  activities  of  23  series  of  drugs  upon  21 
biological  systems,  .“tore  than  2000  biomedical  QSAR  have  been  developed 
involving  more  than  20,000  compounds  and  these  are  stored  in  Hansch 's 
data  bank.  The  strength  of  this  approach  lies  in  the  versatility  cf  these 
substituent  constants  which  apply  to  all  series  of  drugs.  The  importance 
of  the  ideal  partition  coefficient,  Pq,  and  of  the  intercept  in  equations 
in  the  design  of  new  drugs  has  been  demonstrated.  This  method  is  rapid 
but  in  many  cases  only  one  series  of  drugs  can  be  considered  at  one  time. 
Furthermore  wastage  by  metabolism  and  elimination  are  not  accommodated 
by  these  equations, al though  equations  involving  the  same  substituent 
constants  have  been  developed  separately  for  these  processes.  No 
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combination  of  electronic  and  hy*i t ’pnobic  substituent  constants  has  yet  been 
found  to  adequately  predict  tf-  extent  of  hydrogen  bonding  and  its  effect  upon 
the  activity  of  drugs. 

The  method  developed  at  ORES  is  based  upon  evaluation  of  the  three 
dominant  factors  (1.  rate  of  penetration  to  the  site,  2.  rate  of  reaction  at 
the  site  and  3.  rate  of  wastage)  governing  the  degrees  of  biological  activities 
of  drugs  by  in  vitro  physico-chemical  properties  derived  from  suitably  chosen 
model  experiments.  By  applying  the  same  drugs  to  a number  of  different  biologi- 
cal systems  the  significance  of  the  coefficient  of  the  partition  coefficient 
term  in  these  equations  became  apparent.  Mathematical  manipulation  of  equations 
derived  for  the  same  combination  of  drug  families  on  two  different  organisms 
permits  the  calculation  of  the  biological  activities  of  these  drugs  on  one 
organism  from  the  observed  activities  upon  another  organism.  These  equations 
include  a term  to  accommodate  wastage, and  the  experimental  determination  of 
in  vitro  physico-chemical  properties  obviates  the  necessity  of  calculating 
the  effects  of  hydrogen  bonding.  This  method,  involving  the  determination  of 
in  vitro  physico-chemical  properties,  has  advantages  but  the  method  is  more 
tiroe-consumi ng . 

The  Free-Wilsan  method  involves  the  development  of  substituent 
constants  in  biological  activity  units  from  a large  series  of  drugs  in  one 
family  for  one  specific  biological  process.  Addition  of  the  cogent  sub- 
stituent constants  to  the  experimentally  determined  biological  activity  of 
the  parent  drug  gives  the  biological  activity  of  the  derivatives.  These  sub- 
stituent constants  are  all  inclusive, incorporating  into  one  substituent 
constant  the  effects  of  metabolism,  elimination  and  hydrogen  bonding  etc.  The 
disadvantage  of  this  method  is  that  a new  series  of  substituent  constants 
must  be  developed  for  each  series  of  drugs  applied  to  each  biological  process, 
so  it  is  very  time-consuming  and  lacks  the  versatility  of  the  first  method. 

In  all  some  29/  equations  and  207  references  are  presented  in  this 
review  to  illustrate  the  ways  that  these  methods  have  been  used  to  examine  the 
mechanisms  of  biological  processes  and  how  this  information  can  be  used  in  the 
systematic  design  of  more  effective  drugs. 
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SYMBOLS  AND  ABBREVIATIONS 

The  symbols  used  by  Holmes  for  partition  coefficients  in  the 
systems  cyclohexane-water  and  in  1-octanol -water  are  the  reverse  of  those 
used  by  Hansen.  The  same  is  true  for  n and  n'.  For  consistency  in  this 
article  the  following  symbols  will  be  used. 


Partition 

Symbols  used 

Symbols  used 

Symbols  used 

System 

in  this  Paper 

by  Kansch 

by  Holmes 

1 -Octanol -Water 

P 

n 

P1 

P 

;t 

P' 

r.  » 

Cyclohexane-Water 

it 

It 

P 

n' 

- 

n 

A is  the  agonist  in  moles  per  litre  (84  page  351)  cr  per  kg  of  animal 
weight  (84  page  1405). 

Aanai  is  the  analgesic  activity  of  the  agonist  in  moles  per  kg  of  animal 
weight  (84  page  1405). 

An  is  the  stimulatory  activity  of  the  agonist  in  moles  per  litre  in  the 
blepharospasm  test  (84  page  1444). 

Agen  depr  the  9enera^  depressant  activity  of  the  agonist  in  moles  per 
kg  of  animal  weight  (84  page  1405). 

Aresp  depr  the  resPiratory  depressant  activity  of  the  agonist  in  moles 
per  kg  of  animal  (84  page  1405). 

AT  is  the  threshold  stimulatory  activity  in  moles  per  litre  on  the  frog 
flexor  reflex  (84  page  351). 

8R  is  biological  response  usually  recorded  as 

C is  concentration  in  moles  per  litre. 

Ca  is  cats. 

D is  dummy  or  indicator  parameter  (60). 

Ei  is  the  polarographic  half-wave  potential  against  a saturated  calomel 
^ electrode  (84  page  1210). 

E. c.  is  Eactherioiiia  ooLi. 

F. p,  is  a homolytic  equivalent  of  a for  free  radical  reactions  (110). 

E-  is  the  Taft  steric  factor  (27,44). 

I is  the  indicator  or  dummy  parameter  (60). 

I G| q is  the  50J  inhibition  of  growth  caused  by  the  compound  expressed  in 
moles  per  litre  after  17  hours  incubation.  (84  oaae  1344). 

k<-H  is  the  second  order  rate  constant  for  the  addition  of  n-butanethiol  to 
conjugated  heteroenoid  comoounds  (84  page  1105). 


ky  is  the  pseudo  first  order  rate  constant  for  the  hydrolysis  of  conjugated 
heteroenoid  compounds  by  a reverse  aldol  process  (84  page  1104). 

k,,  is  the  rate  constant  for  some  reaction  with  an  aromatic  system  bearing  a 
‘ substituent,  X,  on  the  benzene  ring  and  for  which  k^  is  the  rate  constant 
for  that  reaction  on  the  parent  compound  of  that  family. 

iCu  is  the  equilibrium  constant  for  the  reversible  addition  of  n-C^HySH  to  the 
conjugated  heteroenoid  compounds  (84  page  213). 

Log  P is  the  logarithm  of  the  partition  coefficient  of  the  solute  in  the 
system  1 -octanoi -water  (4).  In  Holmes'  work  (84)  this  is  log  P ’ : see 

the  first  item. 

Log  P1  is  the  logarithm  of  the  partition  coefficient  of  the  solute  in  the 
system  cycl ohexane-water . In  holmes'  work  (84)  this  is  log  P:  see  the 

first  item. 

Log  P"  is  the  logarithm  of  the  partition  coefficient  of  the  solute  in  the 
system  cycl ohexanol -water  (84  page  776). 

Log  P'“  is  the  logarithm  of  the  partition  coefficient  of  the  solute  in  the 
system  ether-water  (84  pages  777  and  781). 

Log  Pq  is  the  ideal  or  maximum  value  of  log  P for  an  agonist  to  cause  maximum 
biological  response  (41).  This  is  derived  from  equations  involving 
(log  P)2  and  log  P by  taking  the  partial  derivative  of  log  BR  with  respect 
to  log  P and  setting  toe  partial  derivative  equal  to  zero  and  solving  for 
log  P.  This  value  is  log  Pq. 

M is  mice. 

MR  is  substituent  molar  refraction  (92). 

u is  dipole  moment. 

n is  the  number  of  compounds  considered  in  a multiple  regression  analysis. 

:i  is  the  increment  (in  log  units)  in  log  P due  to  a substituent  in  the 

system  1 -octanoi -water  (4)(see  n'  reference  84  page  795). 

.T0  is  the  substituent  equivalent  of  log  Pq  (3). 

:: 1 is  the  increment  (in  log  units)  in  log  P’  due  to  a substituent  in  the 
system  cyclohexane-water  (84  page  787  see  n in  Holmes'  nomenclature). 

n ' . „ is  the  increment  (in  log  units)  in  log  P'  in  the  system 

I riier ac  ui  on 

cyclohexane-water  due  to  interaction  between  contiguous  groups. 

(84  page  787) . 

r is  the  correlation  coefficient  for  the  results  derived  from  equations 
developed  by  a multiple  regression  analysis  of  biological  and  physico- 
chemical data. 

R is  rats. 

P.a  is  rabbits. 

rho,  o,  is  the  reaction  constant  in  Hammett's  equation  log  - log  = ,,a. 

It  is  a constant  for  all  substituents  and  depends  only  on  the  reaction 
series  (90,91). 
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s Is  the  standard  deviation. 

S.a . is  Staphylococcus  aioreua. 

S.alb.  is  Staphylococcus  albuB. 

a is  the  substituent  constant  in  Hammett's  equation  log  = oa 

It  is  determined  by  the  nature  of  the  substituent  and  indo jerdent  of 
the  reaction,  the  constant  of  which  is  involved  in  the  equation 
(.90,91  ,92 ,1 03) . * 

c+  is  Brown's  c constant  (105). 

o'  is  Hammett  constant  for  substituted  phenols  and  amines,  etc.  (44,103) 
o*  is  Taft's  o constant  for  aliphatic  compounds  (106,107). 
o ■ is  the  substituent  constant  for  homolytic  or  free  radical  reactions 


TT 

St 


(116,117). 


is  Swain  and  Lupton's  inductive  component 
is  Swain  and  Lupton's  resonance  component 


in  Hammett's  c (108). 
in  Hammett's  a (108). 
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INTRODUCTION 

Search  for  a new  drug  can  be  approached  In  two  ways. 

This  can  be  done  by  the  random  synthesis  and  biological  testing  of  new 
compounds  or  it  can  be  approached  systematically  using  the  knowledge  that 
has  accunulated  from  structure-activity  studies.  Spinks  (85)*  has  examined  the 
problem  of  finding  new  drugs  and  from  his  analysis  he  estimates  that  one  new 
drug  arises  out  of  each  200,000  compounds.  He  also  indicates  that  one  can 
expect  to  find  an  anticancer  drug  cur  of  each  400,000,000  randomly  tested 
compounds.  This  figure  is  about  100  times  the  number  of  known  organic  com- 
pounds. On  the  other  hand,  knowledge  of  the  factors  governing  the  relation- 
ship between  structure  and  biological  activity  is  growing  by  leaps  and 
bounds  (71,72).  Hansen  (70)  has  more  than  2000  biomedical  QSAR  on  more  than 
20,000  organic  compounds  in  his  data  b?.nk  alone.  The  group  at  DRES  applied 
the  same  900  compounds  to  a large  number  of  biological  systems  and  as  a result 
over  500  QSAR  on  some  900  compounds  are  similarly  stored.  In  a search  for  a 
new  drug  then,  one  is  faced  with  the  choice  of  hiring  a very  large  staff  of 
synthetic  organic  chemists  for  a long  period  of  time,  to  say  nothing  of  the 
large  biological  tes+  team  necessary,  or  staking  one's  chances  upon  a QSAR 
approach  or  upon  an  intelligent  combination  of  both. 


References  1-84  cover  quite  completely  the  QSAR  work  of  Hansch,  Holmes,  Lien, 
Fujito  and  thair  collaborators. 
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The  embryonic  concept  which  eventually  bloomed  into  the 
structure-activity  relationship  dates  back  to  1870  when  Crum-Brown  and 
Fraser  (86)  sensed  that  biological  response,  BR,  to  drugs  was  related 
to  their  chemical  structures,  CS.  This  relationship  may  be  expressed 
mathematically  by  equation  1, 

BR  * f{CS) 1 

They  stated  that  it  should  be  possible  to  develop  a calculus  of  SAR 
by  making  small  changes  in  chemical  structure  and  relating  these  to 
BR.  The  obstacle  that  prevented  them  from  realizing  their  dream  was 
the  problem  of  defining  significant  parameters  of  chemical  structure 
in  numerical  terms. 

Two  developments  in  the  first  third  of  this  century  laid 
the  foundation  for  the  development  of  numerical  values  for  evaiuating 
"chemical  structure".  At  the  turn  of  the  century  Heyer  (87)  and 
Overton  (88)  used  oil-water  partition  coefficients  to  quantitatively 
evaluate  the  penetration  of  simple  organic  compounds  which  act  as 
anaesthetics.  It  was  argued  that  compounds  with  high  oil -water 
partition  coefficients  should  be  good  anaesthetics.  While  this  general- 
ization was  true  for  a number  of  compounds,  it  was  not  universally  true. 
Other  factors  must  also  be  operat-ng  in  the  control  of  anaesthetic  activity 
and  these  were  not  apparent  to  Meyer  and  Overton  so  they  were  unable  to 
formulate  an  equation  relating  biological  response,  BR,  to  various  factors 
evaluating  chemical  structure,  CS.  A third  of  a century  of  qualitative 
study  of  the  electronic  effects  of  substituents  by  the  lapworth-Robinson 
School  culminated  in  i 93 5 in  Hammett's  formulation  of  numerical  constants, 
o,  for  electronic  effects  of  substituents  (89,90).  This  simpTe  but 
extremely  important  idea  opened  the  floodgate  for  the  proliferation  of 
specialized  substituent  constants,  j*a+,  a"  and  a-  and  this  led  tc  the 
resolution  of  o into  Inductive  and  resonance  components  (91).  Steric 
constants  were  then  developed  by  Taft  (91).  Molecular  refraction 
(92)  also  provides,  besides  a measure  of  polarizability,  a measure 
of  bulb  volume.  Fol lowing  the  lead  of  Hammett  in  developing  free  energy 
related  subs cituent  constant*,  Hansch  (4)  in  the  early  1960's  developed 


8 


a substituent  constant  n which  is  the  increment  (in  log  units)  in  the  logarithm 
of  the  partition  coefficient,  P,  due  to  a substituent,  X.  This  may  be  expres- 
sed as  in  equation  2.  Log  P and  n adequately  evaluated  the  "random  walk/'  of 

log  Px  - log  ?H  = 2 

drugs  to  the  site  of  action  in  biological  systems  and  also  the  hydrophobic 
binding  of  organic  compounds  to  enzymes,  serums  and  mitochondrial  proteins 
(6,20).  This  permitted  the  expansion  of  equation  1 to  3. 

aBR  = f ( A hydrophobic  + A electronic  + A steric  + A polarizability). 3 

From  extra-thermodynamic  considerations  Hansch  developed  the  general  equations 
4 and  5 relating  biological  response  to  substituent  constants  ar.d  molar 
refraction  The  kj,  k2  etc.  of  equation  4 

log  BP.  = -kj  (log  P)2  + k2  log  P ♦ k3o  + ki»Ec;  + ksMR  + ke 4 

iog  3R  = -fn jH 2 ► m2r:  + m3c  + m4t£  + m5HR  + ^ 

a'vd  the  rrij,  m2  etc.  of  equatJon  5 are  constants.  Under  specified  conditions 
(3)  these  equations  may  he  reduced  to  ones  with  fewer  terms.  In  the  extra- 
thermodynamic  development  of  equations  4 ana  5 it  was  tacitly  assumed  that  the 
rate  of  wastage  (metabolism  and  elimination)  of  the  drug  in  the  biological 
test  system  is  either  zero  or  a constant  value  for  the  family  of  drugs  under 
consideration.  The  word  family  is  used  because  incorporation  of  Hammett's 
sigma,  c,  constant  limits  the  equation  to  the  examination  of  one  famsly  of 
drugs  at  a time.  Equations  have  ueen  set  up  by  Hansch  (18)  relating  1)  the 
biological  response  BR  to  a given  set  of  parameters  and  2.1  tne  rate  of  wastage 
to  the  same  set  of  parameters, which  provides  an  insight  intc  which  member 
of  the  family  provides  the  happy  compromise  between  tne  degree  of  biological 
response  and  the  rate  of  wastage. 

The  advent  of  the  computer  enabled  Hansch,  Holmes  and  others  to 
examine  the  action  of  drugs  on  many  biological  systems  and  to  determine  the 
dominant  factors  governing  the  degree  of  the  biological  response.  Developing 
equations  by  regression  analysis  for  one,  two,  three,  etc.  terms  on  the  right 
■hand  side  of  equations  4 cr  5 gave  correlation  coefficient,  r,  for  the  value 
of  BR,  calculated  from  the  equation  ccinpa"ed  to  the  observed  valuers  we:l 
as  f-e  standard  deviation,  s.  The  enuat^on  with  the  largest  correlation 
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coefficient  and  the  smallest  standard  deviation  would  be  the  equation 
best  fitting  the  data,  The  terns  on  the  right  hand  side  of  the  equation 
reveal  the  dominant  factors  governing  the  degree  of  the  biological 
response.  If  the  coefficient  of  the  term  is  zero,  then  stcric 
factors  do  not  play  a role  in  this  reaction.  If  log  BR  for  the  antifungal 
activity  of  ani 1 ine-phenol  mixtures  against  Candida  albicans  is  measured 
as  log  - where  C is  the  molar  concentration  necessary  to  produce  the 

biological  action,  then  log  ^ is  given  by  equation  6 (35).  Obviously 
the  rate  of  penetration  (log  P)  and  the  rate  of  reaction  (o)  at  the  site 

log  0.555  log  P + 2.193  o - 1.322 

r=  0.982,  s=0.111  6 

of  action  are  the  dominant  factors  governing  ’og  i in  equation  6. 

Holmes  (84),  in  his  approach,  introduced  a term  for  the  rate 
of  wastage  in  the  biological  process  but  has  not  yet  extended  his 
equation  tc  include  terms  for  steric  factors  or  polarizability.  To 
evaluate  the  rates  of  reaction  at  the  site  of  action  he  used  in  vitro 
rate  constants,  kSH,  for  the  addition  of  n-C9H9SH  instead  of  substituent 
constants.  For  the  bacterial  growth  inhibitory  activities,  IG50,  of  the 
conjugated  heteroenoid  compounds,  I(A  = C0CH3,  CQ2C2H5,  C0NH2;  B = 

COCH3,  C02C2Hs,  C0NH2),  the  in  vivo  rates  of  wastage  were  evaluated  by 


I 


\ 


A* 


B 


★ 


In  this  generalized  formula  A will  always  be  considered  to  be  t vans 
to  the  phenyl  qroup  and  B cie  to  it. 


10 


the  'Ik  oi  t.ro  'rates  of  hydrolysis,  ky,  of  the  I compounds  by-  reverse  aldol  pro- 
cess. Multiple  linear  regression  analysis  of  the  data  for  52  compounds  gave 
equation  7. 

log  IQ50  * -0.24  log  P‘*  -0.55  log  k$H  + 1.07  log  ky  + 0.72  . . . 7 

This  50b  inhibition  Of  growth,  IG^,  of  Ctaphyloaocaus  aureus  is  for  17  hours 
at  37°C.  Determination  of  the  stimulatory  activities,  AT  (in  moles  per 
litre),  against  the  frog  flexor  reflex  was  complete  in  5 minutes  so  wastage 
was  not  a significant  factor  and  the  coefficient  of  log  ky  In  equation  8 is 
zero.  The  correlation  coefficient  for  equation  7 was  good  so  Tog  IG55  (calc) 

log  Ay**  = -0.21  log  P'  -0.41  log  -4.59  7a 

from  equation  7 is  about  equal  to  log  IGjq  (obs).  Subtraction  of  equation  7 
from  equation  7a  and  replacement  of  log  IG^q  (calc)  by  log  IG^  (obs)  gives 
equation  7b.  This  and  similar  equations  permits  the  calculation  of  the 

log  Ay  - log  IG5C  (obs)  + G.03  log  P'  + 0.14  log  k^ 

-1  .07  log  ky  -5.31  ....  7b 

biological  activities  of  a family  of  drugs  on  one  organism  frcm  the  observed 
activities  on  another  organism. 

An  alternative  method  for  relating  the  degrees  of  biological 
activity  of  a family  of  drugs  to  de  novo  substituent  constants  was  introduced 
by  Bruice  et  '■.»  (03)  and  more  fully  developed  by  Free  and  Wilson  (94).  Free 
and  Wilson  have  developed  a set  of  substituent  constants  for  a series  of  ten 
tetracyclines  involving  different  substituents  in  three  positions  cn  th? 
parent  ring  system.  By  setting  up  a series  of  simultaneous  equations,  one  for 


* To  avoid  confusion  between  the  terminology  of  Hansch  and  Holmes,  the 
following  symbols  have  been  adopted  throughout  this  review.  P is  the 
partition  coefficient  in  1 -octanol -water , P'  is  the  partition  coef- 
ficient in  the  system  cyclohexane-water  and  P1"  is  used  for  the  system 
ether-water.  This  means  that  P and  P'  are  reversed  in  the  article  by 
Holmes  (84). 

■*  log  i of  equation  6 equals  log  1 - log  C.  This  is  equal  to  -log  C. 

u 1 
Log  At  then  will  be  comnarable  to  -loq  — . 

1 ‘ U 
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each  substituent  in  each  position,  it  was  possible  to  assign  substituent 
constants  (in  biological  units)  so  that  the  biological  activity  is  the 
sum  of  these  biological  substituent  constants  plus  a constant.  This  method 
lumps  into  one  constant  hydrophobic  effects  of  the  substituent  along  with 
electronic,  steric  factors  and  group  interactions.  As  well  metabolism 
and  elimination  are  accommodated  in  these  constants.  The  disadvantage 
of  this  approach  is  that  constants  derived  for  one  set  of  drugs  will  not 
be  useful  for  another  set  causing  a different  biological  response.  This 
is  the  advantage  of  the  free  energy  related  substituent  constants  developed 
by  Hammett,  Taft  and  Hansch  etc.  They  have  been  successfully  applied  to 
over  2000  biomedical  QSAR. 

Strategy  in  drug  modification  is  becoming  a subject  in  its 
own  right  (9S).  Techniques,  independent  of  the  computer,  have  been 
developed  by  Topliss  (41,96,97)  and  Craig  (98)  for  arriving  at  the 
derivative  in  a congeneric  series  of  drugs  with  the  optimum  properties, 
by  examination  of  a minimum  number  of  members  of  the  series.  Cluster 
analysis  (99)  is  being  considered  as  a means  of  minimizing  the  problem  - 
of  collinearity  between  terms  in  multiterm  equations.  As  our  under- 
standing of  drug  action  at  the  molecular  level  advances  through  the 
development  of  biochemistry,  enzymology  (100-102)  and  molecular  biology 
it  becomes  more  urgent  to  systematize  the  interactions  of  organic  com- 
pounds with  as  many  biochemical  systems  as  possible.  The  explosion  of 
data  pouring  from  the  current  journals  is  a waste  of  effort  unless  it  can 
be  retrieved  and  correlated  with  other  data.  This  is  one  of  the  pressing 
problems  in  the  field  of  Structure-Activity  at  the  moment.  The  method 
of  data  management  being  developed  by  Hansch  is  reported  in  three  of 
his  recent  articles  (70-72). 

These  various  approaches  will  be  discussed  in  detail  as  well  as 
the  method  for  deriving  the  ideal  log  P value,  leg  PQ , for  the  member  of 
a series  of  drugs  when  applied  to  some  biological  system.  This  will  be 
used  in  the  systematic  approach  to  the  design  of  more  effective  drugs. 


FREE  ENERGY  RELATED  SUBSTITUENT  CONSTANTS  (92)  ..  . 

Substituent  constants  c and  have  been  developed  for  homogeneous 
reactions  but ,surpri singly, experience  has shown  that  these  apply  equally  as 
well  to  heterogeneous  reactions  of  biological  systems,  for  example -Ej. 
adequately  accommodates  the  steric  interactions  between  substrate  and  - 
biological  system  (40).  Log  P or  n fulfill  a non-specific  role  a, td  a 
specific  role  (37).  By  its  non-specific  role  it  evaluates  the  /'random 
walk"  of  the  drug  to  the  site  of  action,  while  the  specific  role  evaluates 
the  binding  of  organic  compounds  to  the  enzyme  (13).  These  constants  will 
be  discussed  in  general  terms  for  those  not  mathematically  oriented. 

Electronic  Substituent  Constants  (80,91,92,103) 

Determination  cf  rates  of  reaction  between  organic  compounds  or 
substrates  and  the  biophase  at  the  site  of  action  in  a biological  system 
would  be  difficult  or  impossible.  However  the  product  of  the  Hammett  re- 
action constant  p and  the  substituent  constant  a has  enabled  -Hansel)  .to 
obviate  this  difficulty  in  his  approach  to  the  Structure-Activity  problem 
and  he  has  exploited  it  to  the  full. 

To  dispel  any  confusion  as  to  the  meaning  of  electronic  as 
employed  in  this  context  it  will  be  considered  to  mean  the  following.  When 
a drug  reacts  with  a receptor  site  with  the  -formation  of  a covalent  bond, 
the  reaction  will  be  very  sensitive  to  quite  small  changes  in  the  electron 
density  on  and  around  the  atom  involved  in  the  reaction.  Such  electronic 
differences  could  easily  be  small  enough  to  have  little  or  no  effect  upon 
other  parameters,  such  as  partition  coefficients,  but  still  have  a large 
effect  upon  the  covalent  bond  forming  reaction.  Craig  (98)  has  examined  the 
interdependence  of  the  variables  used  in  structure-activi ty  relationships  and 
Hansch  (44)  has  shown  in  equations  8 and  9 that  there  is  a slight  degree  of 
covariance  between  n and  a.  Equation  8 relates  n values  for  meta  substituents 
to  for  similarly  located  groups  while  equation  9 does  the  same  for  para 
substituted  groups. 

n = -1 .84  o +0.70 
m 

34,  r = 0.387,  s = 1.079 


.n  c 


8 


13 


ii  1 -0.89  op  + 0.48 

n 3 37 , r * 0.300,  s * 1 .110 9 

The  Hamnett  reaction  constant  o and  the  substituent 
constant  a are  based  upon  the  observation  of  Hammett  (89)  that  for 
a family  of  benzene  derivatives  there  is  a linear  relationship  between 
the  logarithms  of  the  rate  constants,  k^,  or  equilibrium  constants 
and  the  logarithms  of  the  ionization  constants,  of  the  similarly 
substituted  benzoic  acids.  This  may  be  expressed  mathematically  as 
in  equation  tq  where  the  slope  n is  a measure  of  the  sensitivity  of  the 


log  k^  = p log  + C 10 

logarithm  of  ky  to  electronic  effects  of  the  substituent,  X.  is 
the  ionization  constant  of  tho  similarly  substituted  benzoic  acid. 
When  X = H is  substituted  into  equation  10,  as  in  equation  11,  this 
gives  the  relationship  for  the  parent  compounds  of  these  two  series. 


log  kH  = p log  + C . 11 

Subtracting  equation  11  from  equation  10  gives  equation  12. 

log  i~  ■ p log  . 12 

*H  'Si 

Log  -j-r-  was  defined  by  Hammett  as  a 

'Si 

Kx 

log  ~ = u 13 

Si 

and  substituting  o for  this  term  in  equation  12  gives  equations  14  and  15 
kX 

log  — = log  kx  - log  = pa 14 

kH 


log  k^  = pa  + log  k^ 15 

The  reaction  constant,  p,  by  the  nature  of  the  linear  relationship 
is  a constant  for  all  substituents  and  depends  only  upon  the  reaction  - 


f 

j 


u 

series  being  examined.  The  substituent  constant,  a (in  log  units),  is  by 
definition,  determined  by  the  nature  of  the  substituent  and  is  independent 
of  the  reaction  with  which  k is  associated.  Tables,  of  subs.tituent  constants, 
o, 'appear  in  references  90,  91  and  92,  Substituent  constants,  o,  for  di- 
substltuted  derivatives  can  be  approximated  quite  well  frcm  the  algebraic 
sum  of  the  c constants  for  each  of  the  substituent  groups.  Thus  the  o con- 
stant for  3 ,4-dimethoxy  derivatives  is  the  algebraic  sum  of  for.  the 
methoxyl  group  + op  for  the  methoxyl  group.  While  equation  15  gives  amazingly 
good  correlations  for  a wide  variety  of  reactions  involving  me.ta-  and  para- 

«3.4-d1-0CH,  * "°-V15  * (-°'266)  * -°-153  16 

substituted  benzene  derivatives  (see  ref.  84  pages  215,  220-224,228,236, 

243  and  302),  it  has  its  limitations.  In  the  case  of  ortte-substi tuted 
benzene  derivatives,  steric  factors  are  superimposed  upon  electronic  effects 
so  a different  set  of. constants  (104)  must  be  used.  Equation  15  does  not  hold 
for  aliphatic  compounds. 

There  arc  two  types  of  benzene  derivatives  for  which  Hammett  a 
constants  fail  to  provide  reliable  results.  The  first  one  is  where  a sub- 
stituent containing  lone  electron  pairs  is  coupled  through  a series  of  benzene 
double  bonds  (?)  to  an  electrorv-wi thdrawing  or  accepting  group  (e.g.  CHO,  CM, 

N02  etc.)  as  in  II  and  III,  then  o~  constants  (44)  should  be  used.  A second 
situation  where  a constants  fail  is  when  an  electron-donating  substituent  is  so 
oriented  as  to  reduce  the  positive  charge  at  a centre  on  the  benzene  ring  invol- 
ved in  electrophilic  substitution  as  in  IV  -*•  V. 
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CH3 

IV  V 

Under  these  circumstances,  o+  constants  developed  by  Brown  (105)  and  listed 
in  references  44  and  105  must  be  used. 

Taft  et  al  (106,107)  attempted  to  resolve  a values  into  their 
inductive  and  resonance  components  Oj  and while  Swain  and  Lupton  ( 1 08) 
developed  the  analogous  constants  * and  ft.  Tables  of  these  are  listed  in  ref.  9 

Unlike  the  planer  aromatic  derivatives  with  meta  and  para  sub- 
stituents plots  of  log  k values  for  aliphatic  reactions  against  log  K for 
ionization  constants  of  aliphatic  acids  usually  are  net  linear  (see 
equation  11).  Steric  effects  and  field  effects  in  the  latter  case  are 
not  sensibly  constant  giving  rise  to  a scatter  of  points  in  the  last  plot. 

Taft  (91),  from  a determination  of  the  rates  of  base  and  acid  catalyzed 
hydrolysis  of  esters  of  aliphatic  acids  of  the  type  of  X - CH2C02C2H5 

developed  a o*  constant  and  a steric  constant,  Ec*.  Tables  of  o*  values 
are  listed  in  references  44  and  91,  while  E$  values  are  to  be  found  in 
references  44  and  27.  The  greater  the  steric  effect  of  a group  the 
larger  E<.  is  in  a negative  sense  (66,68). 

Two  sets  of  analogous  constants  o’  (21,109,116,117)  and  ER  (110) 
have  been  developed  for  homolytic  or  free  radical  reactions. 


+ H~ 


★ 


Ej  values  apply  equally  well  to  aliphatic  and  aromatic  substituents  (44). 
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Hydrophobic  Substituent  Constants, n and  n'. 

The  free  energy  related  p and  o constants  of  Hammett  prompted  Hansch 
to  develop  an  analogous  hydrophobic  substituent  constant,  Ji*,  from  partition 
coefficients  in  1-octanol -water,  while  Currie  et  al  (111-114,84)  followed 
with  a set  of  n'*  and  n'.  constants  for  the  systems  cyclohexane- 

water  and  1 -octanol-water . 

Many  partition  syst“ns  have  been  used  (115)  in  the  past  but  recently 
partitioning  has  centred  around  the  1 -octanol-water  (4,84)  and  the  cyclo- 
hexane-water systems  (84).  Comparing  results  for  the  two  systems  reveals  that, 
over  the  central  part  of  the  spectrum  of  solutes  examined,  a linear  relation- 
ship exists  but  divergence  in  opposite  directions  occurs  at  the  two  extremes. 
This  relationship  is  expressed  in  equations  17-20  (84).  Similar  relationships 


log  P - 1 .85  log  P'  -2.53 

n - 51  , r = 0.91 17 

log  P = -0.23  (log  P')2  + 2.98  log  P'  -3.50 

n = 91  , r = 0-93  . . 18 

log  P * -0.088  (log  P')3  + 0.23  (log  P')2  + 2.22  log  P‘  -3.34 

n = 91 , r = 0.93 19 

log  P = -0.016  (log  P’)1*  + 0.045  (log  P')3  -0.115  (log  P;)2 

+ 2.514  log  P'  -3.349 

n - 91,  r - 0.93 20 


between  partition  coefficients  in  the  system  1 -octanol -water  with  cyclo- 
hexanol -water , P"  (84),  and  with  n-butanol -water , P1V  (44)  are  given 


respecti vely 

by  equations  2 

i and 

22. 

log 

P"  s 

: 0.72  log  P + 1 

D.  60 

n » 14, 

r = 0 

.97 

21 

log 

piv 

* 0.70  log  P 

0.38 

n - 57. 

r = 0 

.993,  s = 0.123  ...... 

22 

* To  avoid  confusion,  log  P and  it  refer  to  the  system  1 -octanol -water 
while  log  P'  and  u ' refer  to  the  system  cyclohexane-water.  This  is  just 
the  reverse  of  the  symbols  usad  by  Holmes  in  Structure-Activity  Relation- 
ships (84). 
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It  Is  Icnown  that,,  in  many  cases,  water  penetrates  into  biological 
media  by  pinocytosis*  but  roost  organic  compounds  appear  to  penetrate  into 
biological  media  by  an  iter  ted  process  of  partitioning  between  nonpolar 
and  polar  media.  In  analogy  with  equations  17,  21  and  22  two  partition 
systems  could  be  envisaged  such  as 

System  No.  1 System  No.  2 

octanol  + drug  Biological  protein  + drug 


+ t Kx 


+ f 


water  + drug 


water  + drug 


where  equation  23  would  hold.  Equation  23  relates  the  penetration  of 

log  K2  3 a log  Kj  + b 23 

drugs  into  biological  protein  to  log  P (1 -octanol -water).  It  would  be 
interesting  to  test  this  If  a biological  system  were  available  whose 
response  to  drugs  was  dependent  only  upon  partition  coefficients. 

While  this  is  out  of  the  question,  nevertheless  there  are  many  pure 
enzymes  whose  inhibitions  aredependent  only  upon  these  hydrophobic 
forces  or  partition  coefficients  (44,46). 

By  convention,  the  equilibrium  constant,  P,  for  a solute  in 
the  systan  1 -octanol -water  is  defined  by  equation  24  which  holds  for  the 
process  expressed  in  equation  25. 


[solute  in  1 -octanol]  9i 

P 3 [solute  in  water)  ........ 

Solute  in  water  + Solute  in  1 -octanol  25 

Log  P then  is  a measure  of  the  free  energy  Involved  in  the  reversible 
transfer  of  1 mole  of  solute  from  water  to  1-octanol.  Hansch  (20) 
states  that  "it  has  been  shown  that  the  transfer  of  a hydrocarbon  solute 


* 


Pinocytosis  is  the  penetration  of  compounds  through  holes  in 
cell  membranes  (68). 


from  a nonpolar  environment  to  an  aqueous  one  Is  exothermic  for  aliphatic 
hydrocarbons  and  approximately  athermal  for  aromatics.  The  solubility  of 
these  and  other  organic  compounds  in  water  is  associated  with  the  Targe 
negative  entropy  of  solution  which  is  due  to  the  formation  of  loosely  held 
but  highly  structured  envelopes  of  water  molecules  around  the  apolar 
portions  of  the  organic  molecules  as  they  enter  the  water.  It  is  predominantly 
the  molecular  size  and  shape  which  determines  how  many  molecules  enter  into  the 
structured  sheath  around  the  apolar  portions  of  the  organic  solute  molecules 
and  therefore  determines  the  magnitude  of  the  negative  entropy  of  solution." 
Conversely  the  transfer  of  an  organic  solute  from  water  to  a nonpolar  solvent 
will  strip  the  structured  sheath  of  water  molecules  from  the  solute  molecule 
which  may  result  in  the  generation  of  a weak  (bond  energy  of  about  1 k cal/mole) 
hydrophobic  bond  between  the  nonpolar  portion  of  the  solute  and  the  nonpolar 
solvent.  Hansch  (20)  further  states  that  "the  major  factor  determining  the 
partitioning  of  organic  molecules  between  aqueous  and  organic  phases  is  the 
extent  to  which  they  form  hydrophobic  bonds." 

Hansch  developed  a free  energy  dependent  substituent  constant  H 
as  the  increment  (in  log  units)  in  log  P due  to  a substituent  X when  introduced 
into  the  parent  compound,  H.  Equation  26  expresses  this  relationship.  From  a 
series  of  compounds  values  for  n were  deduced  (92,84)  which,  except  for  minor 

It  = log  Px  - log  PH 26 

variation,  hold  from  series  to  series.  In  fact  log  P for  a compound  can  be 
calculated  from  equation  27. 

log  P = iin  27 

Leo  and  Hansch  (36)  developed  21  equations  by  regression  analysis  relating 
log  P values  for  compounds  in  2:  partitioning  systems  to  those  in  the  system 
1 -octanol -water . When  solutes  that  could  hydrogen  bond  were  combined  in  the 
series  with  those  that  could  not,  the  correlation  coefficient  for  the  equation 
relating  log  P'  to  log  P (analogous  to  equations  21  and  22)  was  poor.  Actually 
two  equations  were  needed,  one  to  accommodate  the  solutes  that  could  hydrogen 
bond  and  another  for  those  that  could  not.  It  is  claimed  that  solvents  such  as 
butanol,  pentanol , cyclohexanol  and  1-octanol*,  which  dissolve  considerable 
amounts  of  water  do  not  require  two  equations. 

* 1 -Octanol  dissolves  2.3  moles  of  water  per  litre  while  1-octanol  is  only 

soluble  in  water  to  the  extent  of  4.5  x 10'3H  (36). 


• ' ■ • (9 

Currie  et  al  (111)  developed  an  analogous  series  of  n*  values 
(see  reference  84  page  787)  for  atoms  or  atom  groups  such  that  log  P' 
for  a compound  is  given  by  equation  28.  For  73  1 compounds  with  B=H 

log  P*  * in'  - 1.30 28 

these  it'  values  gave  log  P*  (calc)  values  that  compared  well  with  the 
observed  values  as  seen  In  equation  29  (84).  Log  P'  values  were  then 


log  P'  (calc)  = 1.00  log  P'  (obs)  + C.04 

n = 73,  r = 1 .00 29 

log  P'  (calc)  = 1.51  log  P'  (obs)  - 1.80 

n = 91 , r * 0.98 30 


calculated  for  91  I compounds  with  two  activating  groups  at  A and  B using 
equation  28, and  the  results  (equation  30)  were  far  from  satisfactory.  The 
correlation  coefficient  was  good  but  the  slope  in  equation  30  is  far  from 
unity.  For  lipophilic  compounds  the  deviation  between  calculated  and 
observed  log  P‘  values  is  as  much  as  + 0.7  log  units  while  that  for  hydro- 
philic compounds  is  as  much  as  -3.7  log  units.  These  deviations  were 
ascribed  to  interaction  between  the  two  contiguous  A and  B groups  in  the 
I compounds.  Accordingly  equation  28  was  expanded  to  31  to  include  a 

. term  (reference  84  page  787)  to  compensate  for  interaction 
between  the  contiguous  groups.  Employing  equation  31  log  P‘  values  were 


log  P'  • in'  ♦ in'1nter4ctl0l) -1-30  31 

calculated  for  325  I compounds  with  one  and  two  activating  groups.  The 
relationship  between  log  P’  (calc)  and  log  P1  (obs)  is  expressed  by 
equation  32. 


log  P1  (calc)  = 0.98  log  P1  (obs)  + 0.08 

n = 325,  r = 0.98 32 


Directing  attention  similarly  to  log  P values  for  the  I 
compounds  in  the  system  1 -octanol -water  it  again  was  found  necessary 
to  introduce  a ninteraction  term  (ref.  84  page  795)  into  equation  27 
to  give  equation  33. 
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W = V + ^interaction  ’ * * / *••'•••  ; 33 

log  P (calc) ’=  1 .07  log  P (obs ) -0.16 

n = 103,  r = 0.97  34 

Log  P values  were  calculated  for  103  compounds  using  equation  33  and  the 
relationship  between  log  P (calc)  and  log  P (obs)  is  given  by  equation  34- 
Benzalacetoacetamide  probably  exists  in  the  hydrogen  bended  forv.  VI. 


Calculation  of  log  P for  this  compound  by  equation  33  (including  n ) 

interaction 

gave  a calculated  value  of  log  P which  differed  from  the  observed  value 
(ref.  84  page  754)  by  0.45  log  units.  (The  deviation  between  log  P‘  (calc) 
and  log  P*  (obs)  was  -1./6  log  units.) 

Folding  of  the  molecule  by  whatever  means  generally  leads  to  a 
more  hydrophilic  compound  than  would  be  predicted  from  either  equation  31  or 
33.  Hansch  and  Anderson  (19)  developed  n constants  in  three  different  ways. 


1 ‘ nX  “ lQgPC6HsCH2CH2CH2X  " log  PC6H5CH2CH2CH3 
2)  nx  3 log  - log  pc6H5CH3 


The  nx  values  developed  by  these  three  methods  all  differed,  however  there 
was  a constant  difference  between  the  values  developed  by  methods  1 and 
3.  These  workers  ascribe  this  to  a folding  of  the  side  chain  in  C6N5ch2LH2Ch2x 
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due  to  interaction  of  the  side  chain  dipole  moment  with  the  n electron? 
of  the  aromatic  ring.  This,  when  reinforced  by  hydrophobic  bonding, 
could  lead  to  a folded  molecule  and  this  leads  to  greater, hydrophi 11  city 
than  would  be  expected.  Currie  at  at  (111)  observed  that  log  P'  for 
vitamin  K.lt  VII,  was  much  lower  than  would  be  predicted  from  Jog  P'for 
2, 3-dimethyl -1 , 4-naphthoquinone.  This  too  may  be  due  to  folding  of  the 
long  side  o.'.tw  due  to  hydrophobic  bonding. 


•ft 


(jH, 


?*> 


VII 


The  phenanthrylene  oxide  bridge  and  the  ethanamine  chain  in 
codeine  Impart  rigidity  to  the  backbone  (rings  A and  B)  of  codeine,  VIII, 
and  neopine,  IX.  The  a7-8  double  bond  in  codeine  leads  to  a puckered 
molecule  while  the  a8"14  double  bond  of  neopine  leads  to  a nearly  planar 
structure  (ref  84  pages  104-109).  Log  P‘“  (ether-water)  for  codeine  is 
+ 0.09  while  that  for  neopine  is  - 0.43.  Obviously  the  folding  of  the 
codeine  molecule  is  blanketing  part  of  the  molecule  from  interaction  with 
the  solvent.  Folding  and  itydrogen  bonding  may  account  for  the  anomalies 


Ns 


CH- 


VIII 


IX 
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present  in  the  log  P values  for  the  hydrochlorides  of  the  benzoroorphans 
listed  in  Table  1 . 

Substituent  Molar  Refraction  Constants,  HR  (92) 

Ec  has  been  successfully  used  in  some  equations  for  evaluating 
steric  effects  in  the  interaction  of  drugs  with  biological  systems,  how- 
ever the  steric  requirements  are  often  of  the  bulk  type  so  Hansen,  for 
want  of  a better  parameter,  has  used  group  molar  refracti vities  which  are 
listed  in  reference  92.  MR,  besides  evaluating  bulk  steric  effects,  is 
also  directly  related  to  polarizability  so  it  should  be  used  with  caution. 
Anotner  limitation  to  the  use  of  MR  in  this  work  is  that  there  is  some 
;)  collinearity  between  MR  and  n (58,67,99)  and  other  variables  (61,70), 

MR  has  been  scaled  down  by  a factor  of  0.1  for  use  by  Hansch  in  references 
v 53,58,59,61,66,67,73  . This  makes  MR  more  equiscalar  w'th  n (53). 


TABLE  1 

LOG  P VALUES  FOR  SALTS  OF  SOME  BENZQMORPHAN 
DERIVATIVES  IN  1 -OCTANOL-WATER 


COMPOUND  LOG  P 


COMPOUND  LOG  P 


-2-Methyl -6, 7-benzomorphan  hydrochloride  -1.65 
2-Methyl -5,9-diethyl -6, 7-benzomorphan  hydrochloride  -0.98 
2-Hydroxy-2-methyl -5-ethyl -6, 7-benzomorphan  hydrochloride  -1.28 
2-Hydroxy-2 ,9-dimethyl -5-ethyl -6 , 7-benzomorphan  hydrochloride  -0.97 
2-Hydroxy-2 ,5-dimethyl -9-ethyl -6 , 7-benzomorphan  hydrochloride  -1.12 
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QSAR  EQUATIONS  DERIVED  FROM  EXTRATHERMODYNAHIC  CONS IDE RATIONS 

Hansch  (1,2,3),  assuming  the  rates  of  metabolism  and  elimination 
to  be  zero  or  constant  and  steric  factors  to  be  insignificant,  has  developed, 
fnom  extrathermodynamic  (69)  considerations,  equation  35  relating  biological 
response  to  hydrophobic  and  electronic  factors.  In  equation  35  biological 
response  lb  usually  recorded  as  log  where  C is  the  molar  concentration 
of  the  drug  eliciting  the  desired  response.  Under  specific  conditions 

log  ~r  = -ki  n2  + k2n  + <3a  + k4 35 


(see  ref.  2 and  3)  this  equation  may  reduce  to  the  sampler  equations 
36-39. 

log  i 3 an  + b 

log  1 = -an2  + bn  + c . . . 

log  * oa  + c 

log  = an  + oo  + c 


36 

37 

38 

39 


If  steric  factors  are  not  insignificant  then  E$  or  MR  can  be  introduced 
into  these  equations  while  a+,  a~  , o*  or  ER  can  be  substituted  for  o 
as  the  situation  demands.  The  data  for  log  £-,n,g,  E<-  etc. can  be  submitted 
to  multiple  regression  analysis  to  derive  the  best  values  for  a,  b,  c etc. 
Deriving  values  for  all  these  equations  and  examination  of  the  correlation 
coefficients  will  indicate  which  equation  best  accounts  for  the  biological 
activity.  Table  2 summarizes  the  types  of  equations  that  have  been 
developed  so  far  while  Table  3 presents  some  of  the  biological  processes 
that  have  been  examined  by  this  method.  This  approach  will  be  illustrated 
by  examining  ten  type  equations  for  QSAR. 

Four  basic  steps  must  be  followed  in  the  development  of  the 
best  equation  to  represent  the  biological  process  and  they  are: 

1)  Equations  are  developed  by  linear  regression  analysis  relating 
the  logarithms  of  the  biological  responses  of  a family  of  drugs 
to  each  of  n,  a etc.,  E$  and  MR.  Even  if  the  correlation  coef- 
ficient and  standard  deviation  are  not  good,  the  correlation 
coefficient,  r,  shows  the  relative  importance  of  P3ch  of  these 
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: TABLE  2 ' . 

SUGARY  OF  EQUATIONS  RELATING  BIOLOGICAL  RESPONSES  TO  VARIOUS 
COMBINATIONS  OF  SUBSTITUENT  CONSTANTS 


Equations  Involving  References  ] Equations  Involving  References 


u (or  log  P) 

7,13,18,23,24, 

25,35,38,40, 

42,45,46,64,68 

n and  y 
a*  and  E<j  . 

62 

15 

0 

34,37 

Es  and  MR 

66 

-f 

0 

21,34,41 

n2,  n and  o 

35 

a~ 

30,40,41 

n2,  a and  o* 
n2 , n and  ER 

21 

28 

0* 

15,61 ,68 

n2,  n and  & 

59 

o • 

21 ,28,41 

(log  P)2,  log  P and 

ES 

68 

ES 

15,22,23,37, 

39,40,51,53 

n2,  n and  0 

40 

n2  and  n (or 
(log  P)2  and  log  P) 

24,35,48, 

50,57,68 

(log  P)2,  log  P and 

PK, 

24 

7.  and  o 

14,15,17,22, 

n , o and  Eg 

15,27 

35,37 

log  P,  o*  and  Ec 

61 

7 and  o 

34,67 

n , £ and  MR 

59 

!!  and  Eg 
a and  MR 

24,37,64 

67 

n,  MR  and  D 

69 

TABLE  3 

SOME  BIOLOGICAL  PROCESSES  SUBJECTED  TO 
QUANTITATIVE  STRUCTURE-ACTIVITY  ANALYSIS 
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Test 

Medium 

Biological 

Action 

Drugs 

No.  of 
Test  Drugs 

Refer- 

ences 

A 

Adenasine  Deam- 
inase (Bovine 
Gut  Mucosa) 

Inhibition  of 

9-(l -alkyl-2' 
hydroxyethyl ) 
adenines 

9 

50 

A&rcbactev 

aerogenes 

Inhibition  of 
Growth 

B-Nitrostyrenes 

14 

84 

Antibody  - 
Antigen  Inter- 
action 

Inhibition  of 

Ben2oate  anions 

22 

26,53 

A.  oleracea 

Kill 

N-Alkylpyridinlum 
Hal  ides 

7 

35 

Quinones 

10 

35 

Arfcacia  Eggs 

Inhibition  of 
Cell  Division 

Barbiturates 

19 

16 

Aspergil Iub 
niger 

Inhibition  of 
Growth 

Benzyl  alcohols 

19 

35 

Isothiocyanates 

13 

35 

B-Nitrostyrenes 

8 

84 

Phenol s 

26 

35 

a,p-Unsaturated 

Ketones 

19 

35 

.ill 

Carboxylate  ions 

8 

35 

Aapeigillue 

Inhibition  of 

Imidazol  ines 

15 

35 

Bolcrni  Growth 
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B 


Bo  oi  Z Xus 
suiitZts 

Inhibition  of 
. Growth 

Rifamycln  B 
amides  and  hydra - 
ztdes 

24  and 

41 

61 

Bacteria 

Inhibition  of 
Luminescence 

Alcohols  and 
Urethanes 

5 and  8 

16 

Bacteria  and 
Fungi  - - 

Antimicrobial 

Activity 

Esters  of  p -Hydroxy- 
benzoic  acid 

4 and  26 

42 

Barnacle  Larvae 

Narcosis 

A1 co ho Is 

14 

44 

Beef  Erythrocyte 
Carbonic  Anbydrase 

Inhibition  of 

4-Subs t -benzene- 
sulfonamides 

12 

50 

Beef  Liver  Mito- 
chondria 

Deamination  of 
Ami nes 

Benzylamines 

7 

24 

Blood  Clot 

Fibrinolytic 

activity 

Salicylic  and 
Benzoic  Acids 

Benzoates,  Sali- 
cylates, anthra- 
nilates  and 
cyclopropane- 
carboxyl  ate  anions 

49 

9 and 
22  - 

33 

64 

Blood  Plasma 

Fibrinolytic 

activity 

Benzoic  Acids 

15 

33 

Salicylic  Acids 

13 

33 

Bctrytia  allii 

Curling  of 
hyphae 

Gri seofulvin 
Analogs 

22 

35 

Inhibition  of 
Growth 

6-Nitrostyrenes 

?0 

84 

Bo  try  tig  oinerea 

Inhibition  of 
Growth 

e-Nitrostyrenes 

8 

84 

Bovine  hemoglobin 

Binding  of 

Miscellaneous 

compounds 

17 

44 

Bovine  Muscle 
Succinate  Oxidase 

Inhibition  of 

Miscellaneous 

compounds 

14 

44 

ZJ 


c 


Candida  albicans 

Inhibition  of 
Growth 

Carboxylate  Ions 

6 ’ 

35 

Diaml nes 

19 

35 

Hydroxybenzoic 

Esters 

7 . 

35 

g-Nitro styrenes 

U 

84 

Pyrimidines 

8 

35,44 

TA3  Carcinoma 
fells  (Mice) 

Antitumor 

Activity 

DiaminopyrimiJines 

10 

68,75 

Carcinoma  of 
Nasopharynx 

Antitumor 

Activity 

N-Acyltri amines 

28 

75 

Cat 

Analgesia 

Morphine  Alkaloids 

i4 

84 

Antagonism  to 
adrenaline 

N-Substituted-2- 
bromoal kylamlnes 

10 

22 

Excitant 

Activity 

Morphine  Alkaloids 

14 

84 

Cattle  Liver 
Amine  Oxidase 

Deamination 
of  Amines 

Primary  Amines 

8 

24 

Cell  Culture 
Nasopharyngeal 
Carcinoma 

Anti  tumor 
Activity 

Acyldiamines 

28 

68,75 

Chloroplasts 

Inhibition  of 
Hill  Reaction 

Anil i des 

25 

14 

Phenyl  dimethyl  - 
ureas 

12 

14 

N-Phenyl isopropyl - 
carbamates 

9 

14 

Cholinesterase 

Inhibition  of 

Diethyl  Phenyl 
Phosphates 

5 

30 

Methyl  carbamates , 

8 and 

15 

Diethyl  phenyl phos-  30 
phates,  Alkyl pbos- 
phonic  acid  esters 
and  Phospho  rami  dates 


28 

- 

; ■ - 

Cholinesterase 
(Human  Plasma) 

Inhibition  of 

RN(CH3)3 

7 

50 

Choi  inesterase 
(Rabbit  Plasma) 

Inhibition  of 

(CH3)3N(CH2)nN(CHj)3 

7 

50 

0 

Di hydrofolate 

Inhibition  of 

Pyrimidines  and 

10  and 

17 

Reductasa 

Triazines 

12 

Triazines 

83  and 
244 

58,73 

Dihydrofolate 
Reductase  {E.aoli^ 

Inhibition  of 

Triazines 

12  and 
15 

17,40 

DihyJrofolate 
Reductase 
(Pigeon  Liver) 

Inhibition  of 

Pyrimidines 

12 

40,50 

01  hydrofolate 

Inhibition  of 

Triazines 

244 

69 

Reductase  (Tucwr) 


E 


Equine  Liver 
Dehydrogenase 

Inhibition  of 

Carboxy amides 

6 

50 

Eaaheriohia  coli 

Inhibition  of 

Chloramphenicols 

10 

28 

Growth 

8-Nitrostyrenes 

23 

84 

F 


Flies 

Synergistic 

action 

1 ,3-Benzdiazoles 
with  Carbaryl 

16 

21 

Frog  Flexor  Reflex 

Stimulation  of 

Benzaldehydes 

13 

84 

Catechol  monomethyl 
ethers 

13 

84 

a-Haloacetophenones 

13 

84 

B'Nitrostyrenes 

29 

84 

Tetrahydro-1 ,4- 

12 

84 

-’l P H + krt/ii j S nAnpc 
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Frog  Flexor  Reflex 

Stimulation  of 

a,e-Unsaturated 

compounds 

46 

84 

Frog  Sciatic  Nerve 

Inhibition  of 

Alcohols 

8 

44 

Fuaarium 

bulbigenwn 

Fungltoxicity 

g-Nitrostyrenes 

8 

84 

G 

2.  cingula  ta 

Inhibition  of 
Growth 

Imidazolines 

15 

35 

Goldfish 

Narcosis 

A1  cohols 

8 

44 

Guanine  Deam- 
inase 

Inhibition  of 

9-Phenylguanines 

32  and 
96 

51,73 

Guinea  Pig 
Complement 

Inhibition  of 

Benzamidines 

108 

53,73 

Benzyl pyridini urn 
ions 

69 

67,73 

Guinea  Pig  Eye 

Blepharospasm 

a,B-Unsaturated 

compounds 

10 

84 

Guinea  Pig  Ileum 

Choi  inerg ic 

E0so 

Choline  Esters 

6 

44 

In  vitro 
Inhibition 

Di phenyl hydrami nes 

30 

27 

Guinea  Pigs 

Anaesthetic 

action 

2-Diethyl  ami  noethyl 
benzoates 

8 

3 

H 

K&uienula  anomala 

In  vivo  Hi s ta - 
mine  Response 

Di phenyl hydrami nes 

22 

27 

Inhibition  of 
Growth 

Phenyl  Methacrylates 

10 

35 

Horse  Heart 
Cytochrome 

Denaturation  of 

Alcohols 

5 

44 

House  flies 

Toxicity  of 

Diethyl  phenyl 
phosphates 

14 

3 

Human  Liver 

Monamine 

N-(pher>oxyethyl ) 

9 

27 

Mitochondria  Oxidase  cyclopropyl  ami nes 

Inhibition 


30 


I 


Influenza  B.  Virus 

Inhibition  of 
Multiplication 

8enziin1da;cles 

15  : 

44 

L 

Letinue  lepideus 

Inhibition  of 
Growth 

Hydroxy  Compounds 

5 

35 

Ammonium  Ions 

6 

35 

Lewis  Lung 
Carcinoma 

Antitumor 

activity 

Nitrosoureas 

8 and 
13 

57,68 

L 1210  Leukemia 
(Mice) 

Antileukemia 

Activity 

Nitrosoureas 

22 

68 

Triazines 

10 

68 

Lobster  Axon 

Resting  Potenti 
Change 

al  Alcohols 

5 

44 

M 

MaavoBporwn 
aarainae  foime 

Inhibition  of 
Growth 

Imidazolines 

15 

35 

M.  au.reus 

Inhibition  of 
Growth 

Rifamycin  B amides 
and  hydrazides 

24  and 
41 

61 

y,  fruetioc  la. 

Inhibition  of 
Growth 

N- Alkyl  ethyl  ene- 
thioureas 

5 

35 

" " (Spores  of) 

Lethal  Dose 

Benzoqui nones 

10 

44 

M.  tuberc'u loazs 

Inhibition  of 
Growth 

Phenols 

14 

44 

Mice 

Anaesthetic 

Activity 

Aliphatic  Ethers 

26 

46a 

Gaseous  anaesthetics  30 

62 

Analgesia 

Morphine  Al kaloi ds 

9 

84 

fintannn1< m of 

Adrenaline 

N-Substitut  ?- 
bromoal  kyl  a 

9 

22 

Antileukemic 

Activity 

Nitrosoureas  and 
Imidazole- 
carbosyamides 

10  and 
22 

41 

Mice 

Convulsant 

Action 

Morphine 

Derivatives 

14 

84 

Lethal  Dose 

Morphine 

Derivatives 

6 and 
14 

84 

Hypnosis  of 

Acetylenic  Alcohols 

8 

44 

Alkyl -aryl -ureas 

23 

44 

Mice,  Rabbits 
and  Rats 

Hypnosis  of 

Barbiturates 

102 

18 

Mon-barbiturates 

73 

18 

Thiobarbiturates 

25 

18 

Monamine  Oxidase 
(Beef  Liver 
Mitochondria) 

Inhibition  of 

9-A1 kyl -8-carbolines 

8 

50 

Isonicotinic  acid 
hydr azides 

6 

50 

it  (1 

(Human  Liver 
Mitochondria) 

Inhibition  of 

A1 cohols 

5 

50 

U ll 

(Mouse  Liver 
Mitochondria) 

Inhibition  of 

Isonicotinic  Acid 
Hydrazides 

6 

50 

<1  II 

(Rat  Brain) 

Inhibition  of 

HOCH2CH(NH2)CONHNHR 

*7 

t 

50 

!4onilia  fructi- 
aota 

Kill 

Qui  nones 

10 

35 

Mosquito  larvae 

Toxicity  to 

Benzoic  Acids 

14 

3 

Mouse  Eye 

Bl epharospasm 

a.B-Unsaturated 

Compounds 

13 

84 

Mouse  TA  3 
Card  noma 

Antitianor 

Activity 

Di ami  nopyrimidines 

9 

75 

N. 

NADH  Oxidation 

Inhibition  of 

Barb i turates 

6 

16 

p 
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Pcnioillium 

ayclopim 

Inhibition  of 
Growth 

Isothiocyanates 

6 

35 

Penicillium  Lumber 
Mould 

Inhibition  of 
Growth 

s-Nitrostyrenes 

8 

84 

Pepsin  (Porcine) 

Inhibition  of 

Alcohols 

8 

50 

Phytophthora  inf cetane 

Kill 

N-Al kyl pyridini um 
Halides 

8 

35 

Plants 

Elongation  of 

Phenoxyacetic  acids 

20  and 
35 

1,2 

Peeudomonae 

aaruginoea 

Inhibition  of 
Growth 

S-Nitrostyrenes 

17  and 
19 

84 

R_ 

Rabbit  Kidney 

Inhibition  of 
Oxidation  of 
Indophenol 

Alcohols  and  Ketones 

8 

44 

Rabbit  Kidney 
Reductase 

Reduction  of 
Acetophenones 

4-Substi tuted 
Acetophenones 

10 

34 

Rabbit  Liver 
Amine  Oxidase 

Deamination  of 
amines 

Primary  and  Secondary 
amines 

10  and 
30 

24 

Rabbit  Liver 

Inhibition  of 
Guanine  Deamin- 
ase 

9- (X-Phenyl )guani ne 

32 

51 

Rabbits 

Glucuronamide 

formation 

A1 iphatic  Alcohols 

11 

24 

Hippuric  Acid 
formation 

Benzoic  Acids 

8 

24 

Progestational 

Activity 

A6-6-Substituted- 

Progesteror.es 

13 

59 

Respiratory 

Depression 

Morphine  Derivatives 

9 

84 

Ost  Drain 

|\U  V Wl 

Inhibition  of 
Oxygen  Con- 

Barbiturates 

10 

16 

sumption 
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Rat  Liver 

Carcinogenic 
activity  of 

Dimethyl aminoa^o- 
benzenes 

41 

3 

Monamine  Oxidase  N-(phenoxyethy1 ) 
Inhibition  cyclopropylamines 

IB 

V 

Rat  Microsomes 

Demethylatlcn 

of 

Tertiary  amines 

18 

10 

Rats 

Antagonism  of 
Adrenaline 

M-Substitut.ed-2- 
bromoalkyl  amines 

n ' 

22 

Antagonism  of 
Adrenaline  & 
Horadrenal i ne 

N,N-Dimethyl-2- 
broraophenethyl  amines 

22 

22 

General 

Depression 

Morphine  Derivatives 

14 

84 

Rat  Tumor 

Inhibition  of 
Guanine  Deamin 
ase 

9-(X-Phenyl )guanine 

18 

51 

Red  Blood  Cells 

Hemolysis  of 

Aliphatic  Acids 

9 

44 

Rodents 

Thyroxine-1  ike 
Activity 

Thyroxin  analogues 

9 

3 

S 

Sheep  Liver  Inhibition  of 

Succinate  Oxidase 

Miscellaneous 

Compounds 

14 

44 

Siaphylccoccua 

albua 

Inhibition  of 
Growth 

S-Nitrostyrenes 

15 

84 

a,B-Unsaturated 

Compounds 

26 

84 

Staphy lococcu8 

aureua 

Inhibition  of 
Growth 

Alcohol s 

9 

4* 

Benzyl  ammonium 
Chlorides 

45 

44 

Chloromycetin 

analogs 

11 

2 

B-Nitrostyrenes 

35 

84 

Penici 1 1 i ns 

8 

9 

Phenoxypeni ci 1 1 ins 

21 

44 

34 

Staphylococcus 

aureus 

inhibition  of 
Growth 

a , e-Unsa tut ated 
Compounds 

35  and 
52 

84 

Staphylococcus 

typhoea 

Inhibition  of  " 
Growth 

Aromatic  Amines 

15 

44 

Phenols 

35  ' 

3 

Streptococcus 
fueea lie 

Inhibition  of 
Growth 

Rifamycin  B amides 
and  hydrazides 

24  and 
41 

51 

S t rep  to  coccus 
hsmohjtiaus 

Inhibition  of 
Growth ; 

Rifamycin  B amides 
and  hydr«2ides 

24  and 
41 

61 

T 

Tadpoles 

Narcosis  of 

Alcohols 

8 

44 

Trichomonas 

vaginalis 

Inhibition  of 
Growth 

S-Nitrostyrenes 

14  ' 

84 

Trichophyton 

asteroides 

Inhibition  of 
Growth 

g-Nitrostyrenes 

8 

84 

Trichophyton 

gypseum 

Inhibition  of 
Growth 

Phenols 

11 

35 

Trichophyton 

intevdigitale 

Inhibition  of 
Growth 

A1 kylpyrazoles 

6 

35 

Carboxylate  Ions 

14 

35 

Diamines 

22 

35 

Tongue 

Rela  Vive  Sweet- 
ness 

m-Nitroanilines 

9 

34 

Trypsin  (Beef 
Pancreas ) 

Inhibition  of 

H2N(CH2)nNH2 

6 

50 

Tyrosine  Hydroxyl- 

Inhibition  of 

Carboxyamides 

6 

50 

ase  (Beef  Adrenal 
Medullae) 


U 


Urease  (Sword  Bean)  Inhi b-' Hr>n  of 


R-CONHOH 


n 


50 


Valyl-t-RNA 
Synthetase 
( E.ooli ) 

Inhibition  of 

RNHt  Salts 

10 

50 

Venturi  a 
inaequalis 

Fungitoxicity 

B-Nitrcstyrenes 

8 

84 

Kill 

N-Alkylpyridiniim 

Halides 

8 

35 

W 

Walker  256 
Carcinoma 

Anti  tumor 
Activity 

Aromatic  Nitrogen 
Mustards 

9 

75 

Aziridinyl-l-phos- 
phinyl  carbamates 

10 

68,75 

Aziridines 

8 

68 

Bis-(i-aziridinyl  )* 
phosphinyl 
carbamates 

10 

75 

Whale  Myoglobin 

Oenaturation  of 

Alcohols 

7 

44 

X 

Xanthine  Oxidase 

Inhibition  of 

9-(X-Phenyl ) 

30 

66,7 

guanines 
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single  parameters  to  the  biological  response. 

2)  Equations  are  next  developed  by  multiple  regression  analysis 
relating  the  logarithms  of  the  biological  responses  to  com- 
binations of  two  of  the  parameters  then  three  etc.  and  the 
correlation  coefficients  examined  to  determine  the  one  with  the 
highest  value  and  the  lowest  standard  deviation. 

3)  In  certain  cases  a doom/.  D,  or  indicator,  I,  parameter  may  prove 
hel pful . 

4)  It  is  difficult  to  draw  any  conclusions  from  the  above  data  if 
collinearity  exists  between  two  or  more  of  the  parameters  in  an 
equation  (70)  so  the  various  parameters  must  be  plotted  against 
each  other  to  determine  if  collinearity  exists  between  any  of 
them. 

Ten  examples  are  now  presented  to  illustrate  this  method. 

1)  Equation  Involving  only  n or  Log  P (46) 

Most  of  the  equations  so  far  developed,  in  which  only  n or  log  P 
are  involved, are  related  to  the  binding  cf  organic  compounds  to  proteins 
(13,20)  and  enzymes  (20,191)  and  the  growth  inhibition  of  bacteria  (46). 

The  binding  of  compounds  to  both  bovine  serum  albumin  (6)  and  bovine  hemoglobin 
(13)  is  dependent  solely  upon  log  P;  electronic  and  steric  factors  having  no 
influence.  A family  of  19  phenols  was  used  in  examining  the  binding,  £• 

(C  = the  molar  concentration  to  produce  a 1:1  complex), to  bovine  serum  albumin 
(equation  40),  while  seventeen  phenols,  amines  and  neutral  compounds  were  used 
in  tne  analogous  study  with  bovine  hemoglobin  (equation  41). 


log  ^ = 0.681  log  P + 2.489 

n = 19,  r = 0.962,  s = 0.133 40 

log  = 0.713  log  P + 1 .512 

n = 1 7 , r = 0.950,  s = 0.160 41 

2)  Equation  Involving  only  a,  a*  or  a* 


Few  examples  have  been  encountered  where  some  biological  process  is 
related  solely  to  a or  o-  etc.  The  reason  for  involvement  cf  only  this  para- 
meter appears  to  be  due  to  steric  interaction  between  the  substituents  of  the 
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compound  with  the  enzyme  or  lipoprotein  membranes  (37).  One  example 
involving  only  o is  in  the  enzymatic  hydrolysis  of  phenyl  sulphates 
(equation  42). 

log  t = 0.930  a + 2,522 

™ n * 1Q>  r » 0.931 42 

McMahon  et  al  (118)  studied  the  enzymic  reduction  of  a series  of 
ten  4-substituted- acetophenones  which  can  involve  resonance  between  the 
substituent  and  the  ketone  group  (see  page  15  ).  Hansrh  developed 
equations  43  and  44  from  the  data  of  McMahon  et  al  for  the  rate  of 
reduction,  kg,  by  rabbit  kidney  reductase  (34), 


log  kQ  = 2.042  a + 1 .173 

n = 10,  r - 0.862,  s = C.467 

log  kQ  = 1 .514  a++  1 .480 

n = 10,  r = 0.914,  s * 0.390 


43 


44 


Obviously  equation  44  (involving  o+  ) fits  the  data  better  than  does 
equation  43  (involving  a). 

In  developing  equations  for  the  bacterial  inhibition  of  a 
number  of  bacteria  by  thirty-six  rifamycin  B amides,  Hansch  (61)  found 
that  equation  45  for  S tap  hy  loco  a cue  hemolytic'us  involved  only  c*. 


log  = -0.93  o*  + 7.83 


45 


n = 36,  r = 0.858,  s = 0.276  

3 ) Equations  Involving  only 

The  superb  work  of  Pauling  and  his  collaborators  (119)  in  the 
hapten-anti  body  interaction  has  been  analyzed  by  Kutter  and  Hansch  (26). 
The  18  haptens  are  all  simple  substituted  benzoate  ions.  Neither n nor  a 
appears  to  exert  any  influence,  however,  the  steric  effects,  E^,  and 
Es  of  the  substituents  in  the  ortho,  meta  and  para  positions  do  as  is 
manifest  in  equation  46. 

log  K = 1.199  E°  + 0.368  E?  -0.467  E^  -1.406 

n = 18,  r = 0.981,  s = 0.217  


4)  Equations  Involving  (Log  P)2  and  Log  P (48) 

In  contrast  to  the  inhibition  of  enzymes  which  is  related  only  to 
Log  ? or  n,  the  equations  for  biological  responses  of  drugs  upon  living 
organisms  usually  involve  a(log  P)2  or  n2  term  (see  equations  35  and  37). 
Experience  has  shown  that  highly  hydrophilic  compounds  in  a series  of  drugs 
generally  have  little  if  any  biological  activity  and  as  the  lipophilicity  . .. 
increases  so.  does  the  response.  Similarly  highly  lipophilic  compounds 
in  the  same  family  of  drugs  show  little  if  any  biological  activity  and  as 
the  hydropnil icity  increases  so  does  the  biological  activity.  The  activity 
from  either  end  of  the  scale  does  not  increase  indefinitely  but  reaches  a 
maximum.  Somewhere  between  ? = 0 and  P - > for  a given  family  of  drugs 
acting  in  a given  biological  system  will  be  an  ideal  P called  PQ  where  the 
biological  activity  will  be  a maximum.  A greater  number  of  molecules  of 
the  drug  with  P^  would  reach  the  site  of  action  in  the  test  interval  of  time 
than  would  drugs  having  another  P value.  The  shape  of  the  above  plot  of 
log  BR  against  log  P leads  to  a parabola:  hence  the  (log  P)2  or  r.2  in 

equations  35  and  37. 

Mansch,  Smith,  Engle  and  Wood  (43,68)  studied  the  antileukemic 
activity  of  a family  of  nitrosoureas  against  BDFi  mice  inoculated  intra- 
cerebrally  with  10**  cells  of  L 1210  leukemia  which  led  to  equations  47  and 
48.  Unfortunately  no  o or  E<-  values  were  available  for  some  of  the 


log  1 = -0.282  log  P + 4.639 

n * 22 , r = 0.886,  s = 0.189 47 

log  Jr  =-0.0568  (log  P)2  -0.0689  leg  P + 4.527 

n = 22,  r = 0.922,  s = 3.163 48 


substituents  in  the  nitrosoureas  so  linear  equations  involving  on-)y  a and  e<. 
could  not  be  developed.  However  the  high  correlation  coefficient  of 
equation  47  indicates  that  tne  biological  activity  of  the  nitrosoureas  are 
highly  dependent  upon  log  P.  The  increase  in  the  correlation  coefficient 
uDon  introduction  of  a (log  P)2  term  is  statistically  significant.  The 
negative  coefficient  of  log  P in  equation  47  is  very  interesting.  The  more 
hydrophilic  the  compound  becomes,  the  greater  is  its  active uy.  This, 
however,  is  only  half  the  problem.  Toxicity  must  also  be  kept  to  a minimum. 


39 


LDio  values,  for  the  nitrosoureas  on  mice  were  used  as  approximate  values 

of  toxicity  and  equations  49  and  50  were  developed  by  the  method  of 

least  squares.  Log  PQ  (see  page  33  ) for  toxicity  is  1 log  unit  higher 

than  that  for  potency  of  the  nitrosoureas.  This  indicates  that  one  should 

log  i = 0.210  log  P + 4.232 

n = 28,  r « 0.737,  s - 0.232  ......  49 

log  J-  = -0.0688  (log  P)2  + 0.0593  log  P + 4.066 

n = 28,  r = 0,809,  s =■  0.206  5Q 

be  able  tc  make  more  potent  nitrosoureas  by  reducing  their  lipophilicity 
and  at  the  same  time  lowering  toxicity. 

5 ) Equations  Involving  n and  a etc. 

Miller  and  Hansch  (17)  have  applied  regression  analysis  to  the 
data  of  Baker  and  Shapiro  (120)  on  the  inhibition  of  dihydrofolate 
reductase  by  pyrimidines.  Equation  51  summarizes  this  relationship. 

log  ^ = +0.457  n -5.820  0 -6.951 

n = 16,  r « 0.903 3 s = 0.741  51 

Turning  next  to  the  work  of  Blanksma  and  Hoegen  (121)  on  the 
sweetness  of  some  substituted-m-nitrani lines,  Xa-**-Xb. 


Xa  Xb 
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Hansch  (34)  developed  equations  52  and  53  relating  the  relative  sweetness 
($S)  to  hydrophobic  (n)  and  electronic  (o  or  c+)  constants.;  Considering 

log  RS  = 1.610  r -1.831  c + 1.729 

n £ 9,  r 3 0.936,  s x 0.282  . ...  . . . . . . 52 

log  RS  - 1.434  n -1.C26  o+  + 1.584 

n = 9,  r = 0.972,  s - 0,190  . . .......  . 53 

the  possibility  of  resonance  between  the  X and  the  W02  groups  in  Xa~+-Xb 
it  is  not  surprising  to  find  that  equation  53  involving  0+  (see  page  15) 
gives  the  higher  correlation  coefficient  and  lower  standard  deviation; 

6)  Equations  Involving  fl  and  £c 

Two  cases  wil 1 be  examined  under  this  heading: 

1)  the  fibrinolytic  activities  of  substituted  benzoate  Ions  (64)- 

2)  the  rate  of  metabolic  change,  MR,  of  m-  and  p-substltuted 
benzylamines  in  the  presence  of  beef  liver  mitochondria  (24). 

The  fibrinolytic  activity  of  the  3-  and  4-substi tuted  benzoic  acids  is 
linear  with  respect  to  log  P but  when  ortho-substituted  derivatives  are 
included  a steric  factor  E$  is  required.  This  led  to  equation  54  for  the 
ortho-substituted-benzoic  acids. 

log  -jL  x 0.48  log  P + 0.44  E^  + 1.36 

n = 16,  r = 0.885,  s = 0.210 54 

In  the  second  case  cited  (24)  the  metabolic  change  (MR)  in  beef 
liver  mitochondria  when  treated  with  meta-substituted  benzylamines  is  linear 
with  respect  to  log  P (equation  55)  but  that  of  the  para-substituted  benzylamines 
is  not  (equation  56).  Insertion  of  an  E$  term  into  the  equation  (eq.  57)  for 
the  para-substituted  benzylamines  markedly  improved  the  correlation  coefficient. 
Combining  the  m-  and  p-substi tuted  benzylamines  gave  equation  58. 


log  HR  » 0.452  log  P + 1.767 

n = 7,  r. * 0.954,  s * 0.0S5 . 55 

log  HR  = 0.256  log  P + 1.229 

n » 7,  r * 0.278,  s - 0.455  . . . ....  . 56 

log  MR  » 0.753  log  P + 0.535  E$  ♦ 0.304  . 

n » 7,  r * 0.787,  s « 0.327  ....  . . . . 57 

log  HR  = 0.623  log  P + 0.683  Es  + 0.554 

n = 13,  r * 0.874,  s * 0.293  58 


7 ) Equations  Involving  n2,  n and  a etc. 

Three  biological  processes  will  be  examined  under  this  heading. 

1}  The  antifungal  activities  of  phenols  against  Aspergillus 
niger  (35). 

2)  The  synergistic  activity  of  some  methyl enedioxy  compounds 
for  the  insecticide  carbaryl  (21). 

3)  The  antibacterial  activities  of  some  chloramphenicols 
against  gram-negative  bacteria  (28). 

A variety  of  antifungal  agents  was  examined  by  Hansch  and 
Lien  (35)  against  a large  number  of  species  of  fungi  and  equations 
developed  by  tha  usual  method.  Of  these,  12  had  equations  of  the 
form  of  equation  36,  sever,  gave  equations  of  the  form  of  equation  39 
and  thirty-two  had  equations  of  the  form  of  equation  37.  Equations 
for  the  action  of  phenols  upon  Aspergillus  niger,  phenyl  methacrylates 
upon  Haneenula  anomala  and  RR'NCSS  Na  upon  Botrytis  cinerea  conformed 
to  the  general  equation  35.  The  equations  for  the  antifungal  action 
against  these  L.hree  organisms  are  respectively  equations  59-61. 

log  ^ = -0.190  (log  P)2  + 1 .859  log  P + 0.627  o -0.092 

n = 18,  r - 0.975,  s = 0.160  .... 

log  «=  -0.120  (log  P)2  + 1.234  log  P -0.880  o + 0,878 


59 


61 
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log  £»  >0.282  (log  P)2  -0.207  log  P -1.531  o + 5,063 

n * 9,  r * 0.921,  s * 0.278  .... 

It  Is  well  known  that  many  methyl enedioxy  compounds  are 
synergists  for  the  Insecticide  carbaryl  (1 -naphthyl -N-metbylcarbamate). 
The  methyl enediaxy  compounds,  XI,  with  the  greatest  synergistic  activity 
were  those  containing  nitro  and  methcxy  groups.  In  nucleophilic  and 

V 

X 

X’ 
y 

XI 


electrophilic  substitutions  these  functions  have  ooposite  effects.  However,  in 
certain  homolytic  (free  radical)  substitutions,  nitro  and  methoxy  are 
amongst  the  strongest  promoters  of  reaction  (116  page  57).  This  prompted 
the  use  of  o-  in  the  equations.  Equations  62-65  were  developed  by 
Hansch  (21)  from  the  synergistic  activities,  SR5*,  reported  by  Wilkinson 
(122).  Hennessy  (123)  from  the  old  classical  approach  to  structure-activity 
postulated  that  the  synergistic  mechanism  involved  the  loss  of  one  nydrogen 


log  SR5  * 0.940a-  + 1.963 

n = 13,  r = 0.603,  s * 0.334  62 

log  SR5  - 0.070  n + 0.917  o-  + 2.050 

n = 13,  r * 0.63d,  s ■ 0.338  63 


★ 


SR5  is  the  synergistic  activity  when  the  ratio  of  synergist  to 
insecticide  is  5:1 . 


log  SR5  3 -0.115  n2  + 0.348  n + 2.146 

n *=  13,  r 3 0.500,  s * 0.380 

log  SR5  = -0.195  n2  + 0.670  n + 1.316  c-  + 1.612 

n * 13,  r 3 0.929,  s * 0.171 
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64 

65 


from  the  methyl anedioxy  group  as  a hydride  Ion.  This  prompted  Hansch  (21) 
to  develop  equations  65  and  67  involving  c+  and  oj  In  place  of  o’. 


log  SR5  = -0.128  n2  + 0.032  ::  + 0.945  jj  + 1.851 

n * 13,  r 3 0.852,  s « C.242  .....  66 

log  SR5  3 -0.113  n2  + 0.374  n -0.166  o+  + 2.184 

n = 13,  r 3 0.532,  s 3 0.392  ...  . . 67 


From  a consideration  of  correlation  coefficients  and  standard 
deviations  equation  65  with  its  free  radlca1  mechanism  would  appear  to 
fit  the  facts  better  than  one  involving  the  loss  of  a hydride  ion. 

Homely tic  constants,  £R,  have  been  developed  by  Yamamoto 
and  Otsu  (110)  (see  page  15  ) and  these  were  usea  by  Hansch,  Kutter 
and  Leo  (28)  in  the  development  of  equations  from  the  data  of  Garrett 
et  al  (124)  on  the  inhibition  cf  S.  ccli  by  a series  of  amphenicols, 
XII.  The  exploratory  equation  68,  developed  by  Hansch  at  al  (2), 


demonstrated  that  activity.  A,  depended  upon  both  electronic  and 
hydrophobic  properties.  The  low  correlation  for  equation  68  prompted 
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Hansch  at  al  (28)  to  employ  the  homolytic  substituent  constant.  Eg  in 
equations  69-73.  Examination  of  the  correlation  coefficients  for  the 

log  A - -0.74  Jl2  + 0.36  n + 1 .82  o„  + Q. 62 

n = 10,  r = 0.924,  s 3 Q.555  ....  68 

log  A " 2.744  t'R  + 0.931 

n 3 8,  r = 0.820,  s * 0.243  ....  69 

log  A = 0.145  n + 1 .289 

n 3 8,  r = 0.317,  s = 0.403  ....  70 

log  A = 0.227  n + 3.069  ER  = 0.769 

n = 8,  r = 0.354,  s = 0.140  ....  71 

log  A ■ 0.187  n + 3.419  £R  - 0.235  a + 0.786 

n = 8,  r = 0.970,  s = 0.127  ....  .72 

log  A = -G. 053  n2  + 0.231  n + 2.865  ER  + 0.846 

n = 8,  r = 0.957,  s = 0.151  ....  73 

two  single  parameter  equations  69  and  70  reveals  that  the  electronic  para- 
meter is  much  more  important  than  the  hydrophobic  one.  While  neither 
correlation  is  good  the  linear  combination  of  these  two  parameters  gives 
an  excellent  correlation.  The  addition  of  a n2  term  does  not  improve  the 
correlation  coefficient  significantly  for  the  limited  range  of  compounds 
exami ned . 

8)  Equations  Involving  n,  o and  (15,27) 

Two  biological  processes  will  be  examined  under  this  heading. 

1)  Inhibition  of  cholinesterase  by  carbamates  and  phosphate 
ester  amides  (15). 

2)  Antihistamines  and  monamine  oxidase  inhibitors  (27). 

Hansch  and  Deutsch  (15)  developed  equations  74-83  from  the  data 
of  Metcalf  and  Fu'xuto  (125)  for  tne  inhibition  of  cholinesterase  by  some 
o-tho-,  met a-  and  para-substi tuted-carbamates , XIII.  From  equations  74-79 
a number  of  salient  features  become  apparent.  There  is  a markedly  different 


4-subs ti tuted  carbamates 

log  = 0.742  n + 3.525 

n = 22,  r = 0.768,  s 

log  1-  3 -1.302  o + 4.202 

n 3 23,  r - 0.404,  s 

log  ~ 3 0.714  n - 0.868  a + 3.486 

n = 23,  r - 0.839,  s 

3-substi tuted  carbamates 

log  J-  = 0.876  K + 4.347 

n = 30,  r = 0.773,  s 

log  J-  = -2.052  a + 5.673 

n = 30 , r - 0.511,  5 

log  ^ ■ 0.784  n - 1 .405  a + 4.618 

n = 30,  r - 0.845,  s 

2-substituted  carbamates 

log  £ 3 0.815  Es  + 4.557 

n = 7 , r = 0. 349 , s 

log  = 1.659  r.  + 4.074 

n = 7 , r 3 C. 566 , s 

log  i-  3 -1  .393  a + 4.887 
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log  £ = 2.79S  U + 4.246  a + 3.845  E$  + 2.542 

n = 7,  r - 0.962,. s = 0.494  ....  83 

character  in  the  activity  of  the  rneta  and  para  isomers.  If  hypothetical 
meta  and  para  Isomers  with  n = 0 are  considered  the  difference  in  activity 
of  the  two  is  1.132  log  units  which  is  a difference  of  over  13  fold. 
Equations  74  and  75  as  well  as  7/  and  78  reveal  that  n accounts  for  a 
much  greater  amount  of  the  variance  in  log  ^ than  does  a.  Moreover  the 
positive  coefficients  associated  with  n in  equations  74,76,77  and  79 
indicate  that  the  more  lipophilic  the  compound  the  greater  is  the  activity 
relative  to  that  of  the  parent  compound.  The  negative  coefficients 
associated  with  a in  equations  75,76,78  and  79  indicate  that  electron 
donating  groups  un  the  phenyl  group  increase  activity  relative  to  the 
parent  compound. 

As  might  be  expected  steric  factors  become  dominant  in  the 
2-substi tuted-carbamates  and  the  negative  coefficient  of  E^'  in  equation  83 
reveals  that  the  smaller  the  substituent  the  more  active  is  the  compound. 

Hansch  and  Deutsch  (15)  also  developed  equations  M to  89,  based 
upon  the  data  of  Fukuto  et  al , for  the  inhibition,  K,  of  cholinesterase 
by  a homologous  series  of  alkyl-2, 4, 5-trichlorophenyl-N-alkylphosphoramides, 
XIV. 


XIV 


The  larger  the  steric  effect  of  the  substituent  the  larger  E<-  becomes 
in  a negative  sense  (66,68,27,44). 
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log  K - 2.709  0*  + 4.490 

n = 8,  r = 0.712,  s = 0.816 84 

log  K * -1 .011  n + 6-567 

n = 8,  r = Q.608,  s * 0.922  85 


log  K = 1.119  Es  + 4.541 

n = 8,  r = 0.875,  s = 0.563  ......  86 

log  K - -3.913  o*  + 2.359  E$  + 4.948 

n * 8,  r = 0.939,  s = 0.438  87 

log  K = 0.001  n + 1 .119  + 4.540 

n = 8,  r = 0.875,  s = 0.617 88 

log  K * 0.007  n - 3.913  a*  + 2.353  E$  + 4.937 

n = 8,  r » 0.939,  s = 0.490  ......  89 


Cf  the  single  parameter  equations,  equation  86  gives  the  best  fit 
thus  demonstrating  the  importance  of  steric  effects.  Comparing 
equation  88  with  86  reveals  that  n contributes  little  if  anything  to 
the  goodness  of  fit  whereas  equations  87  ar,u  84  demonstrate  that  the 
reactivity  parameter,  a*,  plays  a significant  role  The  negative 
coefficients  of  c*  ir.  equations  87  and  89  reveal  that  electron 
donating  R and  R1  groups  promote  activity. 

Fuller  et  at  (126)  determined  the  inhibition,  IG5(>  of  two 
types  of  nonamine  oxioases  hy  N-(phenoxyethyl ) cyclopropyl  ami nes . 

There  is  a marked  difference  in  the  activity  of  isomeric  compounds 
when  the  substituent  is  in  the  meta  or  in  the  para  position. 

Fuller  attributed  this  to  steric  factors.  Rutter  and  Hansch  (27) 
evaluated  the  steric  factors  by  E<.  and  developed  equation  90  for  this 
inhibitory  process. 

log  IG50  = 0.198  n + 1.640  o + 0.702  E$  + 4.153 

n = 18,  r = 0.945,  s = 0.330  . . . 

Kutter  and  Hansch  (27)  developed  a similar  series  of  equations 
for  the  antihistamine  activity  of  a series  of  aryl -substi tuted-di phenyl  - 
hydramines,  XV.  Harms  and  Nauta1  s(  1 27 ) work  was  an  in  vitro  study 


R 


,CH0CH2CH2N(CH3)i 


while  that  of  Ensor  et  al  (128)  was  an  in  vivo  study  on. guinea  pigs.  The 
equations  for  the  two  were  quite  independent  of  11  and  a and  involved  (£$)2  2nd 
£,_  terms.  The  equations  derived  from  in  vitro  and  in  vivo  data  were  surpris- 
ingly similar,  so  much  so  that  both  sets  of  data  were  encompassed  by  a single 
equation. 

9)  Equations  Invol ving  ii2 , n and  D (17,  40) 


Miller  and  Hansch  (17,40)  applied  substituent  constants  to  Baker's 
results  on  the  inhibition  of  dihydrofolate  reductase  by  a homologous  series 
of  1 ,3,5-triazines,  XVI.  In  a number  of  cases  electronic  and  steric  factors 


NH2.HC1 
- R 


were  not  available  for  the  large  R substituents  in  XVI  so  a dunmy  parameter 
was  introduced  into  the  equation.  The  dummy  parameter,  D,  is  assigned  the 
arbitrary  value  of  l.GG  when  a phenyl  group  is  attached  directly  to  the 
triazine  ring  and  a value  of  0.00  when  an  alkyl  group  is  attached  directly 
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to  the  triazine  ring.  This  limited  the  study  to  two  variables.,  n and 
0.  Equations  91  to  93  summarize  the  results.  Incorporation  of  D into 


log  jl  * 0.74  n + 4-25 

n * 15,  r = 0.816,  s = 0.791 91 

log  ~ - -0.23  n2  ♦ 2.12  n + 2.72 

n = 15,  r * 0.887,  s « 0.657  92 

log  J-  = -0.28  n2  + 2.21  n + 0.84  D + 2.68 

n = 15,  r 3 0.928)  s =>  C-553  ........  93 


the  equation  gives  a satisfactory  correlation  coefficient  but  the 
standard  deviation  is  high.  This  is  probably  due  to  the  inability  of 
0 to'exactly  evaluate  electronic  and  steric  factors.  The  positive 
coefficient  of  0 indicates  that  a phenyl  group  attached  directly  to  the 
triazine  ring  leads  to  greater  inhibitory  activity  than  does  an  alkyl 
group  of  equal  1 ipophilicity. 

10)  Equations  Involving  ft,  W and  I (69) 

As  Baker's  work  expanded  to  include  the  inhioitory  activities 
of  more  triazines  upon  di hydrofolate  reductase  from  Walker  256  and  L 1210 
leukemia  tumors,  molecular  refractions,  MR,  for  substituent  groups  were 
being  determined  and  Hansch  and  Si  1 i po  (58)  developed  equation  94  for 
eighty-three  XVII  compounds. 


HH2 


XVII 
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log  £-=  -0.13  a|  + 0.89  n3  + 0.15  MR*  +6-62 


n =83,  r * 0.905,  s * 0.328 


94 


The  subscripts  in  equation  94  indicate  the  position  of  the  substituent,  X, 
in  the  phenyl  group  of  structure  XVII.  At  the  termination  of  Baker's 
classical  work,  Silipo  and  Hansch  (69)  developed,  in  a stepwise  fashion,  a 
series  of  ten  equations  involving  n3  and.  H&+  and  six  indicator  parameters,  I. 
For  example  1-1  made  possible  the  merging  into  one  equation  of  the  activities 
of  the  XVII  compounds  against  di hydrofolate  reductase  from  the  two  types  of 
turners.  The  indicator  value  of  1-1  was  set  at  1 for  Walker  enzyme  data  and 
at  0 for  L 1210  enzyme.  The  complete  equation  for  244  compounds  with 
structure  XVII  is  95. 


0.0243  HR* 


log  J-  = -0.118'  Ii*  + 0.680  n3 

+0.230  MR4  + 0.238  1-1  -2.530  1-2 
-1.991  1-3  + 0.877  1-4  + 0.686  1-5 
+0.704  1-6  + 6.489 
n = 244,  r - 0.923,  s = 0.377  . . 
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Summary 

The  versatility  of  this  approach  to  the  study  of  structure-activity 
relationships  for  bio-medical  processes  is  impressive  and  has  led  to  a 
dramatic  step  forward  in  the  understanding  of  many  of  these  processes, 
however  it  has  some  limitations.  When  log  P values  are  calculated  rather 
than  measured  compounds  with  intramolecular  hydrogen  bonding  should  be 
avoided.  Also  equations  involving  the  electronic  substituent  constants 
a,  c+,  o"  and  a*  limit  the  number  of  families  of  compounds  that  can  be 
included  in  the  analysis  to  one.  Equations  35  to  39  were  developed  on  the 
premise  that  wastage  of  the  drug  by  metabolism  or  elimination  was  not  a 
significant  factor. 

While  no  attempt  has  been  made  to  assess  hydrogen  bonding  or 
loss  of  drug  by  metabolism  and  elimination,  it  is  possible  to  develop 
equations  separately  for  the  rate  of  metabolism  and  uT  e lirni nation  as  well 
as  for  the  desired  biological  potency  and  the  effect  of  log  P or  n on  these 
three  factors.  Hansch  et  al  (18)  has  done  just  this  in  the  study 
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of  the  hypnotic  effects  of  a series  of  barbiturates.  Data  were  taken  from 
the  literature  for  the  hypnotic  effect  of  barbiturates  in  a variety  of 
test  animals  and  under  different  experimental  conditions,  however,  the 
results  were,  inmost  cases,  surprisingly  similar.  The  activity  against 
mice  bears  a parabolic  relationship  to  log  P as  shown  in  equation  96. 

From  data  reported  by  Maynert  and  Van  Dyke  (155)  on  the  per  cent  unchanged 


log  l = -0.438  (log  P )2  + 1 .579  log  P + 1 .926 

n = 13,  r = 0.969,  s = 0.058  96 

leg  % excreted  - -1.235  log  P + 2.695 

n = 1C,  r = 0.957,  s = 0.224  97 

leg  % metabolized  (in  vitro)  - 0.511  log  P + 0.313 

n = 4,  r = 0.987,  s = 0.063  98 

log  % metabolized  (in  vivo)  = 0.634  log  P + 0.599 

n = 3,  r = 0.999,  s = 0.C26 99 


barbiturate'  excreted,  Hansch  et  al  (18)  developed  equation  97. 

Similarly  from  the  data  of  Dorfman  and  Goldbaum  (156)  on  the  per  cent 
of  the  barbiturates  metabolized  by  liver  (in  vitro)  and  by  mice  (in  vivc) 
Hansch  et  al  (18)  developed  equations  98  and  99.  In  contrast  to  equation  96 
equations  97  to  99  are  linear  with  respect  to  log  P but  the  slopes  of  these 
lines  are  in  opposite  directions.  The  optimum  log  P (tha\  is  log  Pq)  for 
equation  96  is  1.80  while  the  % excreted  decreases  wit  ncreasing  lipo- 
philicity  (equation  97).  The  % of  the  barbiturates  metabolized  (equations  98 
and  99),  on  the  other  hand,  increases  with  increasing  1 i pophi 1 ici ty . The 
above  conclusions  are  in  accord  with  general  observations.  Hansch  (44) 
states  that  "ir.  general,  the  more  water  soluble  a molecule  is,  the  more 
rapidly  it  is  eliminated  in  the  urine".  Brodie  et  al  (157)have  demonstrated 
the l ikel ihood  of  a direct  relationship  between  the  rate  of  microsomal 
metabolism  and  the  1 'pophilic  character  of  drugs.  Hansch  e;  al  (10,158) 
have  since  shown  quantitatively  that  as  members  of  a family  of  drugs  become 
mO:‘e  lipophilic-,  other  factors  being  equal,  they  a^e  more  rapidly  destroyed 
by  microsomal  metabolism  (48). 


OSAR  £ QUA T I OHS  INVOLVING  EXPERIMENTALLY  DETERMINED 
CONSTANTS  flN  MODEL  SYSTEMS  (84) 

Holmes  et  at  {84}  obviated  some  of  the  above  difficulties  by  using 
physico-chemical  properties  derived  from  in.  vitro  model  systems  to  evaluate 
the  three  dominant  in  vivo  properties  governing  the  degree  of  stimulatory 
activity  of  the  conjugated  heteroenoid  compounds,  I,  on  the  frog  flexor  reflex 
The  rate  of  penetration  of  the  I compounds  to  the  receptor  site  was  evaluated 


A 

B 

A 

B 

COR 

H 

C0CH3 

COCH3 

co2r 

H 

coch.5 

C02( 

:2h5 

no2 

H 

C32C2H5 

co2c2h5 

no2 

R 

co2c2h5 

conh2 

CN 

CN 

conh2 

conh2 

co2r 

CN 

1? 

conh2 

CN 

/ c^~ 

conr2 

CN 

\ 

<-*33 

/ 

I 


.53 

by  the  logarithms  of  their  partition  coefficients  in  the  systems 
cyclohexane-water  (P‘),  1-octanol -water  (P)  and  cyclohexartol -water  (P"), 
while  the  rates  of  reaction  of  the  I compounds  with  seme  component  of  the 
receptor  were  evaluated  by  in  vitro  second  order  ra  e constants,  k, 
for  the  addition  of  nucleophiles  to  the  conjugated  system  nf  the  I compounds. 
The  rate  of  wastage  in  the  bulk  medium  and  of  the  I compounds  on  the  way  to 
the  -eceptor  site  was  evaluated  by  the  in  vitro  pseudo  first  order  rate 
constants  for  the  hyorolysis  of  these  compounds  by  a reverse  aldol  process 
in  Sorensen  pH7  buffer,  kH  Q,  and  ir.  bacterial  growth  medium  II  of  Schmidt 
and  Moyer  (159),  kM.  The  stimulatory  activity,  Ay  (in  moles  per  litre) 
is  the  threshold  concentration  necessary  for  stimulation  of  the  frog  flexor 
reflex.  The  stimulatory  activity,  A, , was  considered  to  be  proportional 
to  the  product  of  some  power  of  the  above  three  physico-chemical  properties. 
This  can  be  expressed  mathematically  by  equation  100  cr  by  its  logarithmic 
form,  equation  101.  Any  bio-catalytic  effect  that  is  present  in  the 

At  = t(Pa  x kb  x xjj  ) 100 

* rv.  J 

i. 

log  Ay  = a log  P + b log  k + c log  k.^Q  + log  t . . . . 101 

in  vivo  reactions  of  the  I compounds  with  the  receptor  site  would 
be  absent  in  the  in  vitro  reactions  of  the  model  compound  (mirroring 
the  group  cr  groups  involved  in  reaction  at  the  receptor  site)  with 
the  unsaturated  compound.  This  must  be  accommodated  in  equations  100 
and  1 01 . 

A three  step  hypothesis  (Fig.  I)  was  advanced  for  the  chemical 
reactions  involved  in  the  stimulation  of  the  frog  flexor  reflex.  Step 
number  1 involves  the  addition  (XVIII  -*•  XX)  of  some  receptor-nucleophile 
or  -nucleophiles,  XIX,  to  the  conjugated  system  of  the  I compounds.  This 
addition  step  just  generates  the  substrate  for  bio-oxidation  or  bio- 
dehydrogenation (XX  - XXIA  * XXIB)  in  step  number  2.  If  the  receptor  site 
is  regenerated, then  hydrolysis  of  XXIB  to  XXII  (step  number  3)  regenerates 
the  receptor-nucleophile,  XIX, and  XXII  which,  in  turn,  may  suffer  furtner 

am  /-i  r-  m -»  1 "l  r\  f **  a rrma  m f*  c-  r\  Y'  rhcl  af  O uil  t h CnrflO  HI  fa  A 1 T ^ UJA  ^ 
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considered  that  the  bio-oxidation  or  the  bio-dehydrogenation  (step  2) 
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Reactions  Envisaged  in  the  Stimulation  of  Frog  Flexor  Reflex. 
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induces  stimulation  or  initiates  the  step  that  triggers  the  chain  of 
reaction  leading  tc  stimulation. 

If  a reducing  group  were  built  right  into  the  stimulatory 
molecule,  then  neither  the  double  bond  in  the  stimulatory  molecule  nor 
the  addition  step  (step  1)  in  Fig.  I would  be  required  for  stimulation. 
Since  it  is  well  known  that  catechol  monomethyl  ethers  are  readily 
oxidized  to  biphenols  (160-164),  then  stimulation  by  these  compounds  ' 
may  stem  from  the  dimeric  bio-oxidation  or  bic -dehydrogenation  of  these 
catechol  monomethyl  ethers  to  biphenols.  Morphine  and  Its  isomers  are 
readily  converted,  in  the  same  way,  to  pseudo  morphine  (see  ref.  84  page  9) 
and  its  isomers  (165,166),  so  by  analogy  it  was  suggested  that  this  might 
provide  an  insight  into  the  chemical  reactions  associated  with  the 
biological  properties  of  morphine  and  related  compounds. 

Introduction  of  the  bio-oxidation  or  bio-dehydrogenation  step 
would  appear  tc  raise  the  number  of  dominant  factors  from  three  to  four 
for  non-specific  receptor  sites.  If  the  bic-oxldant  at  the  receptor 
site,  however,  were  a very  strong  oxidizing  agent  that  would  oxidize 
or  dehydrogenate  even  weak  reducing  groups  in  the  drug-substrate 
conjugate,  then  bio-oxidaticn  would  not  be  a stimulatory  controlling 
step  even  if  it  is  the  step  inducing  stimulation.  Equation  101  was 
applied  to  the  stimulation  of  the  non-specific  receptor  of  the  frog  flexor 
reflex  by  the  I compounds  and  then  to  the  catechol  monomethyl  ethers. 

This  equation  was  then  successfully  applied  to  the  biological  activities 
of  the  morphine  alkaloids,  the  receptors  for  which  are  stereo-  and 
structurally  -specific. 

The  I compounds  were  demonstrated  to  inhibit  the  growth  of 
four  species  of  bacteria  when  incubated  at  37°C  for  17  hours  in  bacterial 
growth  medium  II  of  Schmidt  and  Moyer  (159).  Combining  equations  for  the 
stimulatory  activities  of  the  I compounds  on  the  frog  flexor  reflex  with 
their  growth  inhibitory  activities  against  bacteria  permuted  the 
calculation  of  stimulatory  activities  on  the  frog  flexor  reflex  from  the 
observed  growth  inhibitory  activities  against  bacteria. 
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Equations  for  Stimulation  of  the  Frog  Flexor  Reflex  and  for 
Inhibition "of  Growth  o?'  Mlcmbrganls'ms  ~ ^ " ' ’ 

Hates  of  hydrolysis  of  the  I compounds  with  two  angular  A and  B 
groups  were.  so  slow  In  Sdrensen  pH7  buffer  and  In  2%  ethanol  -98*  water 
that  it  was  considered  that  wastage  in  the  time  frame  of  the  determination 
of  Ay  would  be  negligible  so  "c"  of  equation  101  was  set  at  zero.  Equations 
were  .developed  by  multiple  linear  regression  analysis  relating  log  Ay  to 
log  P,  log  Pl,  log  P",  log  k$H  {the  second  order  rate  constant  for  the 
■in  vitro  addition  of  n-butanethiol  to  the  I compounds)  and  kg^  (the 
second  order  rate  constant  for  the  in  vitro  addition  of  1 ,3-propanedithiol 
to  the  I compounds)-  These  relationships  are  expressed  in  equations  102  to 
105. 


leg  Ay  * *0.21  log 

P*  -0-41  log  k$H  -4.59 

n 3 46,  r = 0.96,  s = 0.19  . . . . 

. . 102 

log  Ay  3 -0.19  log 

P'  -0.50  log  k2$H  - 4.13 

n 3 12,  r 3 0-95,  s=  0.17.  . . . 

. . 103 

log  Ay  - -0.42  log 

P -0.38  log  kSH  -4.00 

n - 27,  r 3 0.96,  s = 0.19  .... 

. . 104 

log  Ay  «*  -0.42  log 

P"  -0.42  log  k$H  -3.97 

n ° 9,  r * 0.99,  s = 0.10  .... 

105 

The  similarity  In  the  coefficient  of  log  P'  in  equations  102  and 
103  and  of  the  analogous  coefficients  in  equations  104  and  105  prompts 
the  suggestion  that  these  figures  are  characteristic  for  the  biological 
medium  through  whicn  the  I compounds  must  pass.  For  convenience  it  might 
be  termed  the  index  of  penetrability. 

Second  order  rate  constants,  k$H,  for  the  addition  of  n-butanethiol 
to  the  I compounds  have  been  determined  as  have  the  pseudo  first  order  rate 
constants,  k , for  the  reverse  reaction  in  equation  106.  Equilibrium 
constants,  K^,  for  this  reaction  were  then  derived  from  equation  107.  A 
satisfactory  equation  108  was  developed  by  multiple  linear  regression  relating 
log  Ay  to  log  P and  log  for  tins  3-  and  4 substituted  I compounds 
(A=BiC0CH3)  but  it  would  not  accommodate  the  2-substituted  derivatives  nor 


would  It  accommodate  diethyl  benzalmslonate,  I (AaB=C02C2H5}. 
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log  Ay  = -0.26  log  P:  -0.43  log  -2.81 

n = 6,  r = 0.94 


108 


The  rate  constants,  krey»for  the  adducts  of  the  I compounds  (A-B^COC^) 
and  of  (A=B=C02C2H5)  reveal  that  k_ey  is  markedly  dependent  upon  steric 
factors  of  the  A and  B groups  and  of  the  substituent  at  C2.  It  was  also 
seen  t fiat. , exclusive  of  the  ortho-  substituted-3-benzal-2,4-pentane- 
diones,  I ( A=*B= C0CH 3 ) , the  rate  constants, krey,  are  sensibly  constant. 

For  these  compounds,  then,  the  equilibrium  constant,  will  oe  a 


function  of  the  forward  rate  constant,  k 
expressed  in  equation  109. 

log  * 0.91  log  k^H  ▼ 3,89 

n = 6,  r = 0.95 


SH’ 


This  relationship  is 


109 


An  additional  complication  was  considered  when  equations 
for  the  stimulation  of  the  frog  flexor  reflex  by  the  8-ni trostyrenes, 
I(A=*N02,  B-H,  CH3,  C2H-),were  being  developed.  The  observation  of 
Wallace  et  al  (167)  that  under  certain  specified  conditions  nitro 
groups  quantitatively  oxidize  thiols  to  disulphides,  necessitated 


be  ' 

the  consideration  of  this  as  a potential  source  of  wastage  of  the  s-nitro- 
styrenes.  This  reductive  wastage  was  evaluated  by  their  polarographic  half 
wave  potentials,  Ei.  Using  the  A component  of  the  rate  constant  for  the 
addition  of  n-CuHgNH2  to  the  8-nitrostyrenes  afforded  equations  110-113. 


log  Ay  -■*  -0.18  log  P‘  + 0.20  log  A -4.88 

n «■  29,  r = 0.74 .;  . . . 110 

log  At  * -0.16  log  P1  +0.21  leg  A -0-23  log  (E^j  -4.94  . 

n =*  29,  r = 0.74 Ill 

lug  Ar  = -0.72  log  P + 0.24  log  A -3.62 

n * 9,  r = 0.83  112 

log  Ay  = -0.76  log  P + 0.07  log  A + 3.74  log  JE,|  -2.52 

n * 9,  r = 0.88  113 


Since  there  is  no  distinct  improvement  in  the  correlation  coefficients 
wnen  a log  |E^|  term  is  introduced,  it  must  be  concluded  that  reductive 
was+age  of  the  s-nitrostyrenes  Is  not  serious  in  the  stimulation  of  the 
froq  flexor  reflex. 

The  three-step  hypothesis  for  stimulation.  Fig.  I,  suggests  that 
a term  involving  rate  constants  is  not  necessary  fo**  a series  of  catechol 
monomethyl  ethers,  XXIV,  with  about  the  same  critical  oxidation  potential 
(168,169,170).  Equations  114  and  115  were  developed  by  regression  analysis 
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log  Ay 
log  Ay 


-0.38  log  P*  -2.81 

n * 15,  r * 0.83  . . . . 

O.Q9(log  P')2  -0.72  log  P*  -2.60 
n e 15,  r « 0.92  .... 
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The  experimental  conditions  for  determining  Ay  values  for  these  weakly 
active  compounds  were  somewhat  different  from  those  used  for  the 
I compounds  which  may  account  for  the  different  coefficient  of  log  P'. 

The  chemical  reactivity  of  the  I compowds  with  angular 
A and  8 groups  places  these  compounds  in  the  central  portion  of  the 
reactivity  spectrum  of  the  conjugated  heteroenoid  compounds.  With 
the  very  sluggish  members  of  this  series  the  reactivity  parameter 
becomes  the  stimulatory  controlling  step  and  the  coefficient  of  log  P 
or  log  P'  approaches  zero.  At  the  other  end  of  the  spectrum,  where 
the  compounds  are  chemically  very  reactive,  the  partition  parameter 
becomes  the  stimulatory  controlling  step  and  the  coefficient  of  the 
rate  parameter  approaches  zero.  For  the  very  reactive  compounds, 
such  as  the  2-benzal -1 ,3-indanediones , the  rate  of  hydrolytic  wastage 
by  a reverse  aldol  process  becomes  an  important  factor.  As  the  rate 
of  cyanide  addition,  kgN,  increases  so  does  the  rate  of  hydrolysis, 
k^giof  the  benzalmalononitri  les , I(A=B=CN),  in  Sorensen  pH7  buffer. 

This  relationship  is  expressed  in  equation  116.  As  a result  the  -0.80  log 
k^N  term  of  equation  117  for  the  stimulatory  activities  of  the  benzal- 
maloncni tri les  may  be  the  algebraic  sum  of  a stimulatory  component  and 


log  kH  0 = 1 .09  log  kCN  -3.03 

n * 11,  r = 0.96 116 

log  Ay  = -0.26  log  P'  -0.60  log  kCN  -4.77 

n = 11 , r = 0.86 117 


a wastage  component  with  the  stimulatory  component  being  the  larger 
of  the  two.  The  reverse  appears  to  be  true  for  2-benzal -1 ,3-indane- 
diones  (ref.  84,  p.  366). 


The  I compounds  inhibit  the  growth  of  the  bacteria  S.  aureus, 
S.albua,  E.'Coli  and  Aerobaater  aercgenee  but  in  this  review  emphasis 
will  be  placed  on  their  inhibitory  activities  against  s.  aureus.  The 
incuhation  at  37°C  was  for  17  or  18  hours  in  bacterial  growth  medium  II 
of  Schmidt  and  Moyer  (159)>and  not  a matter  of  5 minutes  at  25°C  in 
21  ethanpl  -98S  water  as  in  the  determination  of  A-y.  Furthermore,  Yasnikov 
and  uaivoronskaya  {171}  and  Williamson  and  Witten  (172)  have  demonstrated 
that  amino  acids  and  proteins  catalyze  the  hydrolysis  of  diethyl  benz- 
malonate,  I(AcB*CQ2C.2Hs) , and  related  compounds  by  a reverse  aldol  process 
so  that  it  is  not  surprising  to  find  that,  the  rates  of  hydrolysis  of  these 
compounds,  in  bacterial  growth  medium  II  of  Schmidt  and  Moyer  become  signi- 
ficant in  the  equations  for  inhibition  of  growth,  IG|q  and  1G$®  when  the 
times  are  17  and  18  hours.  The  rates  of  hydrolysis  were  evaluated  by  the 
pseudo  first  order  rate  constants,  k^,  for  the  hydrolysis  of  the  I compounds 
in  bacterial  growth  medium  II.  Regression  analysis  led  to  equations  118-120 
for  the  I compounds  with  angular  A and  B groups. 

For  StapkyloGQGCus  aureus 

log  IG50  * -0.24  log  P'  -0.55  log  It-..  + 1.07  log  ku  + 0.72 


52,  r » 0.88 


log  IG^o  3 -0.51  log  P -0.55  log\SH  + 1.06  log  ky  + 1.44 


n = 40,  r = 0.86 


log  IG50  = -0.17  log  P'  -0.56  log  k2sH  + 0.73  log  ky  -0.22 

n = 11 , r = 0.86  


Again  the  coefficients  of  log  P1  in  equations  118  and  120  are  very  similar 
and  furthermore  the  coefficients  in  equations  118  and  119  are  the  same 
except  for  those  of  the  partition  terms.  The  first  two  terms  on  the  right 
hand  side  of  equations  118-120  must  be  associated  with  the  inhibitory 
process  since, as  log  P,  log  P‘,  log  and  log  k^  become  larger  in  a 
positive  sense.log  IG^q  becomes  larger  (more  active)  in  a negative  sense. 
The  third  term  on  the  right  hand  side  cf  these  equations  must  be  associated 
with  wastage  s1nceias  log  ky  becomes  larger  in  a positive  sense.log  IG5 q 
becomes  smaller  (less  active)  in  a negative  sense. 
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Comparing  equation  121  for  the  inhibition  of  growth  IG$o  {3-. alb) 
of  5.  albue  by  the  I compounds  with  two  A and  8 activating  groups  With 
equation  118  for  S.  aureus  reveals  a change  in  the  coefficients  of  all 
three  parameters. 

log  IG507(5.£rZh)  = -0,15  log  P‘  -0.44  log  kSH  + 0.67  log  ky  -0.75 

n * 27,  r <=  0.70  .......  121 

The  e-nitrostyrenes  and  their  6-alkyl  derivatives  were 
examined  to  study  the  following  factors. 

1)  t)o  dipole  moments  of  the  compounds  evaluate  the  rates  of 
penetration  to  the  site  of  action  as  well  as  or  better  than  partition 
coefficients? 

2)  Is  reductive  wastage  of  the  nitro  compound  (see  page  57  ) 
significant  in  the  17  hour  period  of  incubation? 

3)  What  effect  has  time  of  incubation  upon  the  equation  for 
inhibition  of  growth  of  S.  aureus? 

4)  What  effect  has  the  addition  of  1%  albumin  upon  the  equation 
for  inhibition  of  growth  of  5.  aureu3 ? 

Equations  were  developed  from  the  data  for  S.  aureus  using  for  the  rate 
parameter  the  A and  B components  for  the  rate  constant  for  the  addition 
of  n-C4H9NH2to  the  I compounds  (A=N02,  B=H,  CH3,  C2HS)  and  the  rate 
constants,  k5H>  for  the  addition  of  r.-C4H9SH  to  these  compounds.  A 
log  | E, | term  was  then  introduced  into  these  equations  and  its  effect 
upon  tne  correlation  coefficient  examined.  The  log  P1  term  was  then 
replaced  by  the  logarithm  of  the  dipole  moment  u and  then  combined 
with  log  P1  to  determine  the  effect  upon  the  correlation  coefficient. 

This  is  summarized  in  equations  122  to  133. 

log  IG50  = -0.30  log  P1  + 0.07  log  A + C.70  tog  ky  -1.73 

n = 35  r * 0.90 


i 
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log  16 so  3 -0.28  lag  P'  -0.08  log  A + 3.44  log  |^j  .♦  0.64  log  ky  -1.11 

n = 35 > r * Q .91  . . . . ; 123 


log  IG|o  * -0.26  log  P'  + 0.16  log  B + 0.93  log  ktf  *0.77 

n » 34»  r * 0.91  .....  124 

log  IGso  3 -0.25  log  P*  + 0.13  log  B + 1.10  log  |E,|  + 0.93  log  ky  -0.52 

n = 34,  r 3 0.91  .....  125 

log  IGso  = -0.30  log  P‘  -0.26  log  k$H  + 0.87  log  ky  -0.60 

n = 34,  r = 0.91  . . . . . 126 

log  IG50  55  -0.28  log  P*  -0.16  log  kS(j  + 2.15  log  |E,|  + 0.90.  log. ky  -0.17 

n = 34,  r = 0.92  ^27 

log  Kill  ■ -0.27  log  P*  + 3.00  log  |E,|  + 0.77  log  kg  -0.74 

n-  34,  r - 0.91  . . . . . 128 

log  IG50  = 1.16  log  y -0.10  log  A + 5.15  log  |EJ  + 1.01  log  ky  -0.51 

n = 12,  r = 0.93 129 

log  IG^q  = 1.13  log  y -0.02  log  B + 5.30  log  |E^|  + 1 .13  log  ky  -0.01 

n = 12,  r =0.93 130 


log  IGso  = 1.19  log  u + 0.235  loc  k$H  + 5.62  log  |Ei{  + 0.95  log  ky  -1.09 


n = 12,  r *=  0.93 131 

log  IG50  3 1-15  log  y + 5.07  log  |E^|  ♦ 1-16  log  ky  -0.04 

‘n  = 12.  r * 0.93 132 


log  IG50  10  -0.54  log  P*  + 0.21  log  u -0.51  log  k$H  + 2.94  log  |Ejj 

+ 1.17  log  ky  + 2.13 
n * 12,  r = 0.97 133 

The  procedure  for  the  growth  inhibitory  activities,  IG5®.  reported 
by  Schales  and  Graefe  (173)  for  the  6-nitrostyrenes  and  their  e-alkyl  derivatives 
is  but  a time  modification  of  that  used  abo'e.  The  time  of  incubation  of 
the  S.  aureus  was  18  hours  at  37°C.  Unfortunately  data  on  all  the  compounds 
used  above  were  not  available  for  the  derivation  of  equations  134  to  138. 
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log 

isil 

3 -0.23  log  P* 

+ 0.19  log  A + 

0.90  log  kw  -0.71 

n « 17, 

r ■ 0.94  

* * * 

134 

log 

IGsq 

« -0.19  log  P* 

-0.02  log  A + 

9.62  log  |£^|  + 1. 

06  log 

ky  + 2 

n * 17, 

r * 0.97  . 7 . . . 

. . . 

135 

log 

ir18 

* -0.25  log  P' 

-0.19  loo  k$H  • 

+■  0.82  log  ky  -0.56 

n - 17, 

r = 0.94  

. . . 

136 

log 

TC18 

ig50 

* -0.18  log  P‘ 

+ 0.19  log  kSH  + 10.84  log  |E,j 

+ 0.92 

log  ky 

n - 17, 

r =■  0.97 

. . . 

137 

log 

T p 1 8 

IG5Q 

* -0.18  log  P* 

+ 9.43  log  |E 

+ 1 .09  log  kH  + 

2.34 

n = 17, 

r = 0.97  

« • ■ 

138 

The  average  coefficient  of  log  P‘  for  equations  134  to  138 
is  -0.21  which  does  not  differ  much  from  the  average  coefficient  of 
log  P'  -0.28  for  equations  122-128.  In  general,  the  coefficient  of 
log  ky  and  the  proportional ity  constant  are  larger. in  a positive  sense, 
in  equations  134-138  than  they  are  in  equation?  122-128.  These  two  factors 
contribute,  in  large  part,  to  the  decrease  in  growth  inhibitory  activities 
of  these  confounds  when  the  time  of  incubation  is  increased  from  17  to 
18  hours.  Schales  and  Graefe  (173)  found  that  the  addition  of  1%  albumin 
to  the  bacterial  growth  median  II  of  Schmidt  and  Moyer  reduced  the  growth 
inhibitory  activities,  IG^*3.  of  the  e-ni trostyrenes  and  their  6-alkyl 
derivatives  against  S.  aureus.  In  vitro  studies*  revealed  that  addition 
of  IX  albunin  had  no  significant  effect  upon  the  values  for  log  P‘, 
loa  kcu  and  log  |E,j.  This  amount  of  albumin,  however,  did  catalyze 
the  in  ir'trc  rates  of  hydrolysis  of  the  6-ni trostyrenes  slightly,  while 
the  catalytic  effect  upon  the  6-ni tropropenyl benzenes  was  quite  pronounced 
(see  ref.  84  page  1165).  Rate  constants  ,ky+a , for  the  hydrolysis  of  these 
compounds  in  medium  il  fortified  with  IX  albumin  are  reported  in  ref.  84 
page  1165.  Substituting  log  ky+a  for  log  kw  gave  equations  139  to  141  for 
log  I G 5 against  S.  aureus. 


it 


C.E.  Lough.  Unpublished  data. 
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*,  r,n  T f ^ 

log  IG50 


X 


log  IG 


18+a 

so 


log 


ipis+a 

1«50 


-0.16  log  P*  -0.33  log  A + 11.51  log  |E^|  -0.24  log  ky+a 


-2.10 

n = 17*  r = 0.81 139 

-0.12  log  P'  + 0.39  log  k$H  + 11.81  log  l^j  + 0.23  log  ky+a 

-1.24 

n » 17,  r = 0.80 140 

-0.18  log  P'  + 3.82  log  lE^  + 0.21  log  ky+a  -2.43 

n^=  17,  r = 0.70 141 


If  the  albumin  is  adsorbed  on  the  surface  of  the  bacteria  or 
permeates  into  the  interstices  of  the  cell  wall,  then  this  will  alter  the 
permeability  characteristics  of  the  cell  wall.  These  characteristics  will 
trend  towards  those  cf  the  gram  negative  bacteria  such  as  E.  aoli, with  high 
protein  content  in  the  cell  wall.  The  coefficients  of  log  P*  in  equations 
139-141  are  more  positive  than  those  in  equations  134-138  as  are  those  in 
the  equations  for  E.  aoli. 

The  antimicrobial  action  of  the  B-ni trostyrenes  and  their  B-alkyl 
derivatives  against  S.  albus , t.  aoli,  Aerobaater  aeryoger.ee.  Pseudomonas 
acrugi>\08d.  Trichophyton  aeteroidec , Bo  try  ti  3 allii  , E.  oir.erea, 

Venturi  a inaequalis  , Fu  Barium  bulbigenum , Penicillium  Lumber  Mould, 
Aspergillus  r.iger,  Candida  albicans  and  Tr'ichomonas  vaginalis  were  treated 
similarly  in  ref.  84. 

Equations  for  Various  Biological  Activities 
of  the  Morphine  Alkaloids 

a-Isomorphine , XXVI,  B-i somorphi ne , XXVII,  and  y-i somorphi ne , XXVIII, 
are  isomers  of  norphine,  XXV.  Morphine*,  XXV,  and  a-i somorphi ne , are  diastere- 
omers  as  are  e-and  y-i somorphines . 


* For  the  structures  of  the  morphine  alkaloids  see  the  treatises  by 
Manske  and  Holmes  (165)  and  by  Bentley  (177). 


S- Iscmorphi ne ; XXVII.  and  morphine,  XXV,  are  structural  isomers,  as  are 

y- i somorphi r.e  , XXVIII,  and  a-i somorphi ne , XXVI.  Since  the  receptors 

involved  in  the  biological  actions  are  stereo-  and  structural ly- 

specific  (ref  S4  page  460),  a judicious  choice  of  derivatives  most  be 

made  from  the  large  group  (129  compounds)  of  morphine  alkaloids  for 

which  biological  activities  have  been  reported  by  Kreuger,  Eddy  ana 

Sumwalt  (174).  Analgesic,  ''an,ig>  excitant,  Aexcit>  and  emetic,  Aernet> 

activities  have  Deen  determined  against  cats,  Ca*,  respiratory  depressant, 

A . , activities  nave  been  determined  against  rabbits,  Ra , 

resn  dcpr 

while  general  depressant,  A n ,jepr»  activities  have  bean  determined 
against  rats,  R.  As  well  convulsant,  A , activities  and  lethal 
doses,  A.  . , nave  been  determined  on  mice,  M.  More  recently, 
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Drs  L.J.  Sargent  and  A.E.  Jacobson*  have  determined  the  analgesic  activities 
of  2 G morphines  and  codeines  on  mice  by  the  standard  arid  reproducible 
(175)  method  of  Eddy  and  Leisibach  (176).  The  results  from  this  method 
have  the  added  advantage  that  all  the  data  have  been  subjected  to  probit 
analysis. 

On  the  hypothesis  outlined  on  page  53,  these  biological 
activities  should  be  dependent  upon  four  dominant  factors,  1)  the  rate  of 
penetration  to  the  receptor  site,  2)  the  fit  factor  of  the  alkaloids  with 
the  receptor,  3)  the  rate  of  reaction  at  the  receptor  site  and  4)  the  rate  of 
wastage  or  of  transformation  to  another  morphine  derivative  prior  to 
reaching  the  receptor  site. 

Since  structural-  and/or  stereo-specific  receptors  are  involved, 
the  selection  of  morphine  alkaloids  and  derivatives  to  be  Included  in  the 
regression  analysis  is  much  more  difficult  than  for  non-specific  receptor 
sites  and  involves  a numbar  of  arbitrary  decisions  as  to  which  compounds 
have  structures  that  can  be  considered  to  be  sensibly  the  same.  Any  trans- 
formation of  the  morphine  molecule,  such  as  acetylation  of  the  C6 -hydroxyl 
to  6-acetylmorphlne  or  methylation  of  the  C3-hyd!-oxyl  group  to  codeine, 
represents  a change  of  structure.  Catalytic  hydrogenation  of  the  ethencid 
linkage  probably  alters  the  conformation  of  morphine  and  its  derivatives 
which,  in  turn,  may  alter  the  relationship  between  the  rate  cf  penetration 
and  the  logarithm  of  the  partition  coefficient.  These  changes  will  also 
influence  the  fit  factors  of  the  alkaloids  relative  to  the  topography  of 
the  receptor  which,  in  turn,  may  alter  the  rates  of  reaction  at  the 
receptor  site.  As  a result,  many  compromises  must  be  accepted  if  a 
series  of  alkaloids  is  to  be  subjected  to  regression  analyses.  This  and 
the  larger  error  in  the  determination  of  the  biological  activities 
probably  contribute  to  the  poorer  correlation  coefficients  for  the  plot 
of  calculated  activities  against  observed  values. 
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The  five  criteria  used  in  the  selection  of  the  alkaloids  for 
analysis  were  as  follows: 

1)  the  structure  of  the  alkaloid  must  be  firmly  established; 

2)  the  relationship  between  the  rate  of  penetration  to  the 
receptor  site  and  the  logarithm  of  the  partition  coefficient, 
P"‘  (in  ether-water ) .must  be  constant  for  all  members  of  the 
series; 

3)  the  fit  factor  of  the  alkaloids  with  respect  to  the  topography 
of  the  receptor  site  must  remain  sensibly  constant; 

4)  the  mechanism  of  the  chemical  reactions  involved  at  the 
receptor  causing  the  biologies1  response  must  be  the  same 
for  all  members  of  the  series; 

5)  the  rate  of  wastage  of  the  members  of  the  series  must  be 
sensibly  constant. 

These  restrictions  limited  the  study  to  an  analysis  of  simple 
derivatives  of  morphine,  n-isomers  and  their  dihydro  derivatives. 

Equations  first  order  and  second  order  with  respect  to 
log  P1"  were  developed  by  regression  analysis  for  the  biological  tests 
reported  by  Kreuger,  Eddy  and  Sumwalt  and  by  Sargent  and  Jacobson. 

This  analysis  is  summarized  in  equations  142  to  153a. 


,05  Aanal g-M*  * -°'64  loS  P'" '5'67 

n * 9,  r = 0.95 142 

log  Aanal  M = -0.067  (log  P"‘  )2  -0.63  log  P‘"  -5.63 

n = 9,  r = 0.96 143 

108  Aanalg-Ca  = -°'35  109 

r?  = 14,  r = 0.72 144 


* See  list  of  symbols  and  abbreviations. 
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l°g  Aana]g_Ca  = 0.005  (log  P'“  )2  -0.36  log  Plu  -5.99 

n c 14,  r = 0.72  

145 

l09  Aexcit-Ca  " -°-57  log  p'"  '5-96 

n * 14,  r = 0.86  ..... 

146 

log  Agxcit_Ca  = 0.216  (log  P"*  }2  -0.51  log  Pm  -6.04 

n = 14,  r = 0.88  

147 

’“3  A«sp  depr-Ra  = -0-3’  ’°S  P'“  ‘6-95 

n = 9,  r = 0.92  

148 

1o9  Aresp  depr-Ra  = -0-062  (1og  P'H  )2  "°-83  1og  P'"~6-94 

n = 9,  r = 0.92  

149 

loo  A Hon.  n ” -1  ■ 04  log  P1"  -5.10 

3 aen  depr-R 

n = 14,  r = 0.94 

150 

109  Agen  depr-R  = °‘276  (1°9  P‘''  )2  "1'G9  log  P"'  '5-21 

n = 14,  r = 0.95  

151 

109  Aco„v-M  ■ -0-88  ’°9  P'"  -3-19 

n - 6,  r = 0. 96 

15? 

!og  Alc_h  - -0.79  log  P1"  -3.12 

n = 6,  r = 0.95  

153 

!og  Al[)_m  = 0.129  (log  P"'  )2  -0.84  log  P,fc  -3.13 

n = 6 , r = 0 . 96  

153a 

This  would  indicate  that  the  relationship  between  the  logarithm  of  the 
partition  coefficient  and  the  rate  of  penetration  holds  for  an  iterated 
process  involving  transfer  of  a drug  through  many  cells  and  across  a 
multiplicity  of  membranes.  Secondly,  with  judicious  caution,  this  method 
for  examining  drug  action  can  be  applied  to  systems  that  involve  stereo- 
and  structural ly-specific  receptors. 
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Calculation  of  Biological  Activities  on  One  Organism 
from  Observed  Activities  on  Another 

Equations  have  been  developed  relating  the  biological 
activities,  A (in. moles  per  litre),  of  the  agonists, I,  to  several 
of  their  in  vitro  physico-chemical  properties.  These  equations 
may  be  represented  by  the  generalized  equation  154  where  0 is  some 
in  vitro  measure  of  the  agonist  as  an  oxidizing  agent.  Similar 
equations  have  been  developed  relating  the  bacterial  growth  inhibitory 
activities,  IG  (in  moles  per  litre). of  these  agonists  to  the  same 
in  vitro  physico-chemical  properties.  These  equations  can  be 
represented  by  the  generalized  equation  155.  In  these  equations 
the  coefficients  of  some  of  the  terms  may  be  zero.  Subtraction  of 
equation  155  from  154  leads  to  equation  156.  If  the  correlation 
coefficient  for  log  IG  (calc)  vs  log  IG  (obs)  is  good,  then  log  IG  (obs) 
can  be  substituted  for  log  IG  (calc)  in  equation  156  to  give  equation 
157. 


log 

A 

(calc)  = 

a log  P'  + b log  k + c 

log 

0 + “ '°9  kwast  * 

log  t 

. 154 

log 

IG 

(calc) 

* a'  log  P1  + b' 

log  k 

+ c 1 

log  0 + d 1 log  k .+  log  t' 

• • 

155 

log 

A 

(cal c)  = 

log  IG  (calc)  ■+ 

■ U-a') 

log 

P ‘ + (b-b 1 ) log  k 

♦ (c-c‘ ) log  0 

+ (d-d' 

) log 

W * 109  (/t' 

156 

log 

A 

(calc)  = 

log  IG  (obs)  + 

(a-a ' ) 

log  P 

1 + (b-b ' ) log  k 

+ ( c- c 1 ) log  0 

+ (d-d1 

) log 

k„ast  " '°9  t/V 

. 157 

The  correlation  coefficient  for  log  A (calc)  vs  log  A (obs)  for 
equation  157  would  not  be  expected  to  be  as  high  as  that  for  equation 
154  since  two  experimental  errors  are  introduced,  namely  in  the 
determination  of  log  A (obs)  and  log  IG  (obs).  In  some  cases  within  a 
single  series  of  compounds  the  two  experimental  errors  may  deviate  in 
the  same  direction  from  the  true  value,  while  in  others,  the  deviations 
may  be  in  opposite  directions.  A second  factor  that  may  lead  to  reduced 
correlation  coefficients  is  if  one  term  such  as  c log  0 in  equation  154 
is  the  algebraic  sum  of  two  temis  representing  two  biological  processes 
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such  as  stimulation  and  wastage,  while  the  analogous  term  in  equation  155 
represents  but  a single  process.  Collinearity  between  terms  in  the 
equation  will  also  add  complications.  Equation  157  has  been  applied  to  the 
data  for  a number  of  organisms  using  I compounds  with  angular  A and  b 
groups  and  I compounds  where  A = NOj  and  B 3 H,  CH3  and  C2H5.  This  met nod 
applies  equally  well  to  the  calculation  of  biological  activities  of  the 
morphine  alkaloids  from  their  observed  activities  upon  another.  This  is 
illustrated  by  equations  158  to  175  (Fig.  II). 

The  analgesic  activities  of  a few  morphine  and  codeine  derivatives 
on  man  have  been  reported  (175).  Equations  have  been  developed  from  these 
limited  data  and  the  ultimate  goal  of  calculating  the  biological  activities 
of  drugs  on  man  from  their  observed  activities  on  test  animals  appears  now 
to  be  within  reach. 


log  A . n = log  A _ 

a gen.depr.-R  a excit.-Ca 
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Suren  ary 

The  locus  of  points  for  the  activities  of  a series  of  drugs 
against  one  organism,  when  plotted  against  those  of  another  organism, 
are  often  far  from  linear  and  the  reason  for  this  becomes  apparent  from 
a comparison  of  log  Ay  values  for  the  I compounds  with  angular  A and  B 
groups  with  those  for  the  stimulation,  Ag^**(in  moles  per  litre)of  the 
mouse  eye.  The  relation  of  observed  log  AgW  and  log  Ay  values  is  given 
in  equation  176.  To  determine  the  extent,  if  any,  of  wastage  in  the 

log  A&m  = 0.62  log  Ay  -0.72 

n = 13,  r = 0. 51 176 

stimulatory  process  on  the  frog  flexor  reflex  by  the  I compounds  with 
angular  A and  B groups,  the  pseudo  first  order  rate  constant,  ky,  for 
hydrolysis  in  bacterial  growth  medium  and  log  |E||  were  inserted  in 
equation  102  to  give  equation  177.  Equation  178  is  the  analogous 
equation  for  stimulation,  of  the  mouse  eye. 

log  Ay  = -0.20  log  P‘  -0.45  log  k$H  -1.47  log  I | -0. 03  log  ky  -4.61 

n * 22,  r = 0 98  177 

log  A^  = -0.05  log  P'  + 0.05  leg  k$H+1.96  log  j E^!  -0.78  log  ky  -6.74 

n = 13,  r = 0.77 178 

These  equations  are  quite  different.  In  both  cases  the  coefficient  of  the 
log  ky  term  is  negative  so  it  is  associated  with  stimulation*  rather  than 
wastage.  Furthermore  the  coefficient  of  log  P‘  in  equation  177  is  much 
larger  in  a negative  sense  than  that  in  equation  178.  Subtraction  of 
equation  177  from  178  and  replacement  of  log  Ay  (calc)  by  log  Ay  (obs) 
(r=0.98)  yields  equation  179.  The  correlation  coefficient  for  the  plot  of 

log  AgM  = log  Ay  (obs)  + 0.15  log  P‘  + 0.50  log  k$H  + 3.43  log  |EJ, 

-0.75  log  ky  -2.14 

n = 13,  r = 0.92 179 


* Log  ky  in  equation  178  is  a better  model  than  log  for  evaluating 
the  rate  of  reaction  at  the  receptor  site. 

**  These  values  determined  by  sequential  blepharospasm  test  by  Mr.  B.J.  Wenner 


75 


log  Ag^  (calc)  vs  log  (obs)  is  high  (0.98).  Comparing  equation 
179  with  176  it  will  be  seen  that  the  extra  terms  In  equation  179 
account  for  the  improved  correlation  coefficient.  The  larger  coef- 
ficient of  log  P1  in  a negative  sense  in  equation  177  compared  to 
that  in  equation  178  indicates  that  the  frog's  leg  favours  the 
penetration  of  lipophilic  (oil-soluble)  compounds  to  a greater  extent 
than  does  the  mouse  eye.  The  equations  for  the  B-nitrostyrenes  and 
their  6-alkyl  derivatives  against  S.  aureus  and  B.  allii  clearly 
Illustrate  this  point.  Equations  180  to  183  were  developed  from  the 
data  of  McGowan  at  al  (178)  for  log  1G100.  The  average  coefficient 


log  IGl()o  = 0-58  log  P'  -0.02  leg  A -0.82  log  ky  -8.67 

n = 6,  r - 1.00 180 

log  IG ! oo  = 0.6Q  log  P‘  -0.02  log  A + 5.32  log  |E1(I 

-0.35  log  ky  55.56 

n a 19,  r s 1.00 181 

log  IG100  = 0.59  log  P'  -0.15  log  + 6.78  log  |E^| 

+0.03  log  ky  -3.48 

n = 22,  r = 1 .00 182 

log  IG100  = O.GO  log  P:  + 5.32  log  ] E ^ | -5.46 

n = 19,  r = 1 .00 183 


of  log  P'  for  equations  122  to  128  related  to  5.  aureus  is  -0.28  while 
that  for  equation  18C-183  for  B.  allii  is  + 0.59.  The  log  IG10C  value 
for  3-methoxy-4-hydroxy-fi-nitrostyrene  (log  IG100  = -5,39)  against 
8.  allii  reveals  that  it  is  10  times  more  active  than  6-ni trostyrene 
(log  1G100  = -4.38)  even  though  it  is  chemically  much  less  reactive 
than  6-nl trostyrene.  This  difference  in  activity  stems  mainly  from 
the  relatively  greater  hydrophilic  property  of  3-methcxy-4-hydroxy- 

5- ni trostyrene  (log  P*  = + 0.04)  compared  to  that  of  e-ni trostyrene 
(log  P'  = + 1.80).  The  high  lipophilic  property  (log  P*=  +3.45)  of 

6- nitrobutenyl benzene , I(A=N02,  B=C2H5),  contributes  materially  to  the 

low  activity  (log  IG100  = -2.85)  of  this  compound. 

For  those  equations  involving  a hydrolytic  wastage  term  to 
hold  for  the  I compounds,  the  activities  of  the  two  derived  fragments 
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must  be  either  very  low  or  a constant  value.  The  log  IGjq  values  of  the 
3*  and  4-substituted-benzaldehydes  against  S.  aureus  were  one  order  of 
magnitude  lower  than  for  the  I compounds  with  two  angular  A and  B groups 
andthelrlog  IG^q  values  were  sensibly  constant.  The  activities  of  the 
other  hydrolytic  fragments  were  low. 

Methods  have  been  developed  for  calculating  log  P,  log  P* 

(page  19),  log  k^  and  log  ky  purely  from  structure  on  the  blackboard 
which  permits  the  calculation  of  log  AT  and  log  IG50  values  purely  from 
their  chemical  structures. 

Covariance  or  col  linearity  between  the  various  in  vitro  para- 
meters is  a complication  that  makes  interpretation  of  the  mechanism  of  the 
in  vivo  process  difficult.  One  such  covariance  is  expressed  mathematically 
in  equation  116.  So  any  substituent  introduced  into  tiie  phenyl  group  of 
the  I compounds  to  enhance  reactivity  will  also  enhance  the  rate  of  wastage. 
Maximum  stimulatory  activity  of  the  I compounds  with  A=B=CN  or 


will  result,  then,  from  a happy  compromise  between 


enhanced  activity  due  to  increased  rate  of  reaction  on  the  one  hand  and  loss 
of  activity  due  to  an  increased  rate  of  wastage  on  the  other.  This  optimum 
can  be  approached  from  either  of  the  two  extremes  (84  page  393). 

1)  Introduce  stronger  activating  groups  at  A and  B in  the  I compounds 
to  increase  the  rate  of  reaction. 

2)  Introduce  large  steric  groups  into  the  2 and  6 position  of  the 
2-benzal-l  ,3-indanediones  to  slow  down  the  rate  of  reaction. 

The  parameters  used  in  these  equations  have  been  plotted  against  each  other 
to  test  for  collinearity  or  covariance.  These  equations  are  catalogued  in 
Fig.  Ill, 


Eauations  Reflecting  Col  1 ineari ty  Between  Various 


so 


W NOVO  SUBSTITUENT  CONSTANTS 

A completely  empirical  quantitative  approach  to  structure- 
activity  relationships  was  first  presented  in  1956  by  Bruice  et  at 
(93)  ard  later  taken  up  by  Free  and  Wilson  (94),  Ban  and  Fujita  (179), 
Beasley  and  Purcell  (180),  Craig  (181)  and  Clayton  and  Purcell  (182). 

The  basic  assumption  In  this  method  is  that  substituent  effects  (in 
biological  units)  when  added  to  tne  biological  activity  of  the  parent 
compound  give  the  biological  activity  of  the  derivative.  If  the 
biological  activities  of  a series  of  drugs  belonging  to  one  family 
are  available  a series  of  simultaneous  equations  can  be  set  up  with 
one  unknown  for  each  substituent  in  a given  position.  Solution  of  n 
simultaneous  equations  for  n unknowns  provides  substituent  constants 
in  biological  units  for  the  series  of  drugs  in  question  and  for  that 
specific  biological  process.  Craig  used  69  2-pnenylu qui nol i nes  to  yield 
69  equations  for  their  antimalarial  activity.  The  substituent  con- 
stants, so  derived,  are  listed  iri  references  181  and  44.  Hansch, 

Silipo  and  Steller  (65)  applied  this  method  to  the  study  of  the  inhibition 
of  di hydrofclate  reductase  by  105  2,4-diamino-5(3,4-dich1orophenyl)-6- 
substituted  pyrimidines,  XXIX.  A correlation  coefficient  of  0.920  was 
obtained  for  these  105  compounds. 


Cl 


XXIX 
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The  advantage  of  this  approach  is  that  all  properties  of  each 
substituent  (hydrophobia  electronic,  steric  and  other  more  subtle  pro- 
perties) are  encompassed  in  each  substituent  constant.  Modern  chemistry 
has  long  recognized  the  great  importance  of  electronic  and  steric  factors 
of  substituents  on  rate  and  equilibrium  processes.  As  a result  it  would 
seem  highly  unlikely  that  constants  derived  for  one  set  of  drugs  on  one 
biological  process  would  be  applicable  to  another  set  of  drugs  applied 
to  a different  biological  process.  Group  interactions  complicate 
the  situation  still  further. 

DISCUSSION  OF  RESULTS 

Hansch  clearly  demonstrated  that  the  degree  of  biological 
activity  of  drugs  is  dependent  upon  three  dominant  factors: 

1)  hydrophobic  effects, 

2)  electronic  effects. 

3)  steric  effects. 

This  method,  however,  takes  into  account  neither  wastage  by  metabolism 
and/or  elimination  nor  the  effects  of  hydrogen  bonding.  The  method  of 
Holmes  takes  the  rate  of  wastage  and  hydrogen  bonding  into  account  but 
it  has  not  been  extended  yet  to  include  steric  factors.  The  Free- 
Wilson  method  lumps  all  factors  into  one  substituent  constant. 

The  rates  of  penetration  or  the  hydrophobic  effects  are 
most  important  in  determining  the  degree  of  biological  activity  of 
drugs.  For  equation  102,  log  P ' accounts  for  73%  of  the  variance  in 
the  variable  log  Aj.  Hansch  has  shown  that  the  partition  coefficient 
is  associated  with,  and  evaluates  two  effects. 

1)  The  non-specific  effect  which  evaluates  the  random 

walk  through  the  biological  medium  to  the  site  of  action  (41). 

2)  The  specific  effect  which  is  a measure  of  the  weak  hydro- 
phobic  bonding  between  drug  and  enzyme  (48). 

The  non-specific  role  of  the  hydrophobic  effect  is  dramatically  illustrated 
by  the  equation  derived  by  Hansch  et  al  (11)  of  the  quantitative  deter- 
mination by  Solway  (181)  of  the  amount,  C^,  of  mono-  and  di -substituted 
benzeneboronic  acids  reaching  the  brain  of  mice  15  minutes  after  injection. 
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By  regression  analysis  of  the  data.Hansch  ez  al  (11)  showeo  that  for  pene- 
tration into  the  brain  the  best  fit  is  given  by  equation  210. 


log  Cb  = -0.540  n2  + 0.765  n + 1.505 

ti  = 25,  r = 0.915.  s = 0.214 210 

Inclusion  of  c in  the  regression  analysis  had  no  effect.  In  this  case 


evaluates  the  random  walk  of  the  benzeneboror.ic  acids  to  the  brain. 

Many  examples  of  the  specific  hydrophobic  effect  are  known. 

Many  equations  have  been  developed  by  Hansch  relating  log  P or  n to  the 
oinding  of  drugs  to  pure  enzymes.  These  equations  are  generally  *irst  order 
with  respect  to  log  P or  n as  is  the  case  for  binding  molecules  of  42 
organic  compounds  to  bovine  serum  albumin  (44).  Equation  211  relates  the 
concentration,  C,  necessary  for  the  formation  of  a 1:1  complex  to  log  P 

log  ~ = 0.75  log  P + 2.30 

n = 42 , r = 0.960,  s = 0.159 211 

The  series  of  42  compounds  is  not  a homologous  series  but  is  comprised  of 
phenols,  anilines,  aliphatic  alcohols,  ketones  and  naphthalene  as  well  as 
rigid  molecules  such  as  hydroxyadamantane  and  camphorquinone . Molecules 
which  hydrogen  bond  wall,  such  as  phenols  and  alcohols,  are  accommodated 
equally  as  well  by  equation  211  as  are  molecules  with  little  or  no  ability 
to  hydrogen  bond  such  as  naphthalene,  azobenzene  and  chloroni trobenzene  and 
those  prone  to  form  charge  transfer  complexes  all  fit  the  above  equation. 
Hansch  and  Anderson  (19)  presented  evidence  that  hydrophobic  bonding  between 
a long  hydrocarbon  side  chain  and  an  aromatic  ring  contributes  to  maintaining 
such  a molecule  in  a folded  form  thus  materially  altering  its  log  P value. 
Hydrophobic  bonding  is  considered  to  be  an  integral  part  of  the  partitioning 
process  (20).  Blanketing  of  part  of  a molecule  by  large  bulky  groups  or  by 
virtue  of  part  of  a molecule  being  held  in  a rigid  position  can  prevent 
solvation  of  parts  of  the  molecule  thus  leading  to  anomalies.  The  following 
two  examples  clearly  illustrate  this  point.  Blanketing  has  also  been  observed 
by  Currie  et  al  (113,84  pages  96  and  106). 

ni-C3H7'  ::n-C3H7  = '°'13 

_ u -n  r w = -0.22 
X,  D"Cun0 
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Log  Pq  and  iiq 

It  has  been  demonstrated  many  times  that  equations  for  the 
biological  activities  of  drugs  in  living  organisms  (the  benzenebaronic 
acids  in  mice;  equation  210) are  second  order  with  respect  to  log  P, 
while  in  vitro  studies  with  enzymes  often  lead  to  equations  that 
are  linear  with  respect  to  this  parameter.  In  many  instances  this 
has  led  to  failure  in  the  design  of  inhibitors  for  living  animals 
from  the  data  for  optimum  activity  for  in  vitro  inhibition  of  enzymes. 
Baker  found  many  extremely  potent  inhibitors  [in  vitro)  for  dihydro- 
folate reductase  which  were  quite  ineffective  on  leukemia  cells.  The 
locus  of  the  points  in  the  plot  of  a second  order  equation  with  respect 
to  log  P (equation  210)  is  a parabola, while  that  for  a first  order 
equation  with  respect  to  log  P (211)  is  a straight  line.  Equation  211 
mathematically  states  that  binding  power  increases  indefinitely  with 
increase  in  log  P,  while  in  the  case  of  equation  210  an  optimum  concen- 
tration is  reached  at  a certain  n value  and  then  the  concentration 
reaching  the  brain  progressi vely  falls  off  with  increase  in  log  P. 

The  optimum  value  of  n , namely  Hq,  for  equation  210  is  obtained  by 
taking  the  derivative*  of  log  Cb  with  respect  to  n and  setting  this 
equal  to  zero.  Log  F for  the  parent  benzenebcroni c acid  is  about  *1.6 

d log  C. 

2.  =-2  x 0.54  n + 0.765 

dr 

2 x-  0.54  n + 0.765  = 0 

H,  = 0.71 
0 

so  log  P^w^l  be  1.6  + 0.71  - 2.31.  The  explanation  of  Baker's 
observations  is  shown  graphically  in  Fig.  IV.  The  activity  Gf  the 
in  vitro  inhibitor  with  log  P = +3.90  would  be  very  high  (straight 
line)  but  the  same  compound  would  have  a very  low  activity  in  the 


* 


If  c is  included  in  the  equation  then  the  partial  derivative 
with  respect  to  n must  be  taken. 
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in  vivo  experiment  (parabola)*. 


FIG.  IV 


* The  in  vi tro  action  of  antihistamines  on  guinea  pig  ileum  (127)  and  the 
ir.  vivo  action  of  these  compounds  upon  guinea  pigs  (128)  is  the  exception 
to  the  above  statement.  Kutter  and  Hansch  (27)  have  developed  equations 
relating  these  two  activities  to  substituent  constants.  Both  activities 
are  independent  of  r and  ii2  terms  and  are  dependent  solely  upon  E<. . 

Hence  the  rate  of  penetration  to  the  site  of  action  in  the  in  vivo 
experiments  is  not  a controlling  factor  in  determining  tne  biological 
response.  In  fact  both  sets  of  data  correlate  very  well  on  one  equation. 


AT 

OP 


Equations  212  to  227  derived  by  Hansch  et  al  (18)  from  others’ (182-189) 
data  for  the  hypnotic  activity  of  barbiturates  and  non-barbiturates 
permitted  the  calculation  of  log  PQ  values  (Fig.  V).  These  data 
were  obtained  in  a number  of  laboratories  from  1923  to  1949.  Not  only 
was  hypnosis  defined  in  different  ways  such  as  ED  and  MED50  but  some 
workers  ustd  rabbits,  others  dice  and  still  others  rats.  Even  the 
techniques  used  in  the  same  Taboratory  varied  over  the  years.  The 
average  log  Pq  for  equations  212  to  227  is  + 1.98  and  this  does  not 
differ  markedly  from  that  for  the  localization  of  the  benzeneboroni : 
acids  (log  PQ  = 2.31)  in  the  brain  of  mice  (equation  210).  Similar 
results  have  been  obtained  for  the  inhibition  of  growth  of  Gram 
negative  (equations  228-234)  and  Gram  positive  (equations  235  to  243) 
bacteria.  The  constancy  of  log  Pn  values  for  each  of  these  biological 
systems  led  Hansch  at  al  (18)  to  conclude  that  “other  factors  being 
constant,  se^s  of  congeners  (families  of  drugs)  acting  by  the  same 
mechanism  on  Che  same  receptor  sites  should  have  the  same  log  Pq 
values."  So  if  a new  series  of  hypnotics  is  being  developed  then 
tne  member  of  the  series  with  log  P = m +2.0  should  have  the  maximum 
hypnotic  activity.  While  log  Pq  or  nQ  is  a very  useful  function.it 
must  be  emphasized  that  it  is  an  empirically  determined  constant  and 
that  other  factors  may  influence  the  value  of  log  Pq.  For  example 
side  effects  such  as  metabolism  and  elimination,  which  are  log  P- 
dependent,  may  materially  influence  the  value  of  log  Pq.  For  many 
nan-polar  functions  ::  and  molar  volume  (MR)  are  col  linear  so  increase 
in  n may  influence  the  fit  factor  between  drug  and  receptor. 

It  is  surprising  how  many  of  the  equations  in  Fig.  V have 
about  the  same  intercept  (the  average  value  of  the  intercept  is  + 1.685) 
in  view  of  the  different  test  animals  used  and  the  diversity  of 
techniques  employed.  In  certain  instances  intercepts  can  be  very 
helpful  in  the  design  of  drugs. 


Equations  and  Log  PQ  for  Hypnosis  by  Barbiturates  (Eqs  212-219)  and  non- 
Barbiturates  (Eqs  220-227). 


Equations  and  Log  Pq  for  inhibition  of  Growth  of 
Gram  Negative  Bacteria. 
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Slopes  and  Intercepts  In  the  Equations 

Comparison  of  slopes  and  intercepts  for  equations  involving 
different  families  of  compounds  on  the  same  biological  test  system 
or  different  families  of  compounds  on  different  test  systems  provides 
an  insight  into  the  similarity  of  the  biological  processes  and  their 
sensitivity  to  response  by  these  chemicals.  The  slopes  of  equations 
244  to  264  are  a measure  of  the  sensitivity  of  the  biological  systems 
to  the  hydrophobicity  of  the  chemical  series  examined.  Equations 
244-264  have  been  developed  by  Hansch  at  al  (42)  from  data  gleaned 
from  the  literature  on  the  antimicrobial  activities  of  various 
esters  of  p-hydroxybenzoic  acid  against  a number  of  Gram  positive  and 
Gram  negative  organisms  as  well  as  some  fungi.  These  equations  are 
listed  in  Fig.  VIII.  The  average  slope  of  the  Gram  positive  organisms 
is  0.863,  while  those  for  the  6rair.  negative  organisms  and  the  fungi  are 
respectively  0.540  and  0.518-  There  is  a marked  difference  in  the 
sensitivity  of  the  Gram  positive  organisms  to  the  hydrophobic  properties 
of  the  esters  of  p-hydroxybenzoi  c acid  from  those  of  the  other 
two  groups.  Within  the  Gram  positive  orv_ni$rn  there  is  little  constancy 
in  the  intercept. 

The  magnitude  of  the  intercept  is  determined  by  the  sensitivity 
of  the  biological  system  and  the  chemical  reactivity  of  the  active  groups 
in  the  families  of  drugs  under  examination.  If  the  same  family  cf  drugs 
is  applied  to  two  different  biological  systems  then  the  intercept  reflects 
the  sensitivity  of  response  of  these  two  systems  to  this  family  of 
cnemicals.  Selection  of  an  ester  of  p-hydroxybenzoic  acid  with  a P*  1 
or  log  P » 0 and  substitution  of  this  value  in  equations  244  and  248 
yields  equations  265  and  266. 


of  p-Hydroxybenzoic  Acid  agcinst  Gram  Positive,  Gram  Negative 
Organisms  and  Fungi  (42). 


S.  awreuft 
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eubtitie 


log  ~ = -1.011 
log  ^ * +0.1 97 
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These  two  organ! 'ns  have  the  same  sensitivity  to  hydrophobi city , but  the 
sensitivities  (equations  265  and  266)  of  these  two  organisms  to  this  chemical 
are  quite  different.  This  must  be  associated  with  the  stereo-electronic 
character  of  the  group  in  this  ester  that  is  involved  in  the  inhibition  of 
growth  of  these  two  organisms. 


Collinearity  of  Parameters 

Any  significance  attached  to  relationships  discerned  in  multi- 
parameter equations  of  the  form  of  type  equations  4-10  (pages  36  to  50) 
must  be  tempered  by  a consideration  of  possible  covariance  between  parameters*. 
Covariance  between  various  in  vitro  physico-chemical  properties  used  in 
determining  log  Ay  and  log  IG50  values  has  been  discussed  on  pages  76  to  79. 


* Another  factor  which  can  lead  to  false  conclusions  is  if  the 
activities  of  metabolites  are  of  the  same  order  of  magnitude 
as  that  of  the  drug.  Usually  metabolites  are  quite  inactive, 
however,  an  exception  to  this  is  the  primary  metabolites  of 
A8-  and  A9-tetrahydrocannabinol  (192). 
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For  certain  collections  of  substituents  Hansch  has  demonstrated  that  a 
correlation  exists,  to  a degree,  between  r.  and  n (37,41),  n and  W 
(53,60,61  ,67,70),  o*  and  t<-  (37)  and  MR  and  E^  (61).  A correlation 
matrix  was  then  presented  by  Kansch  et  al  (99)  for  eight  parameters 
involving  90  substituents. 

APPLICATIONS  OF  QSAR 

OSAR  has  not  yet  been  developed  and  refined  to  the  stage  where 
it  can  forecast  a "new  find"  or  drug  for  some  specific  biological  p-ocess 
but  does  provide  a systematic  and  direct  route  to  the  development  of 

the  most  effective  drug  in  the  drug  family  of  the  "new  find".  Ex  mination 

of  the  equations  developed  for  an  already  known  series  of  drugs  upon  this 
biologica1  process  reveals  the  relative  importance  of  hydrophobic  electronic 

and  steric  factors  and  permits  the  calculation  of  the  optimum  log  Pq 

value  and  in  many  cases  Es-q  values  for  the  oest  drug  in  the  "new  find". 

The  above  factors  also  shed  some  light  upon  the  topography  of  enzyme 
surfaces.  Combining  this  information  with  a good  knowledge  of  chemistry 
should  lead  to  a series  of  specifications  for  a drug  to  fit  the  require- 
ments for  the  specific  biological  process.  It  would  seem  to  the  author, 
then,  but  a short  step  to  successful  drug  design. 

Predictions  for  More  Active  Congeners 

Log  Pq  and  signs  of  coefficients  of  the  various  terms  if!  the 
generalized  equations  267  and  268  provide  an  insight  into  the  trends 
that  will  lead  to  optimum  activity. 

log  -j~  = -log  C = an  + bo  -i  rEj  + log  t 267 

log  I G 5 q = a'  log  P + b*  log  + c'  log  k^  + log  t . . 268 

If  the  coefficients  "a"  and  "b"  of  equation  267  are  positive  then, 
as  n and  get  larger  in  a positive  sense, log  C will  get  larger 
in  a negative  sense  (more  active).  If  the  coefficient  "c"  of  equation 
267  is  positive  then  as  E<.*  gets  larger  in  a negative  sense  log  C 

t A positive  sign  associated  with  o indicates  that  it  is  an  electron 
-withdrawing  substituent  (Cii,  N02  etc)  while  a negative  o value 
implies  that  it  is  an  electron  donating  substituent  (Ch^.C^O.vO^^N  etc) 

* As  the  substituent  gets  larger  becomes  larger  in  c negative  sense. 
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becomes  larger  in  a positive  sense  (less  active).  So  with  positive  coef- 
ficients in  equation  267,  increases  in  1 ipophillclty  and  stronger  electron 
withdrawing  groups  will  enhance  biological  activity,  while  a decrease  in 
size  of  the  substituent  will  produce  a trend  in  the  same  direction.  The 
signs  of  the  coefficients  in  equation  268  reveal  whether  this  term  is 
associated  with  the  biological  process  or  with  wastage  (see  page  60). 

Moreover  the  larger  "a"'  is  in  a negative  sense  the  more  the  biological 
media  favour  the  penetration  of  lipophilic  compounds. 

If  the  equation  for  the  biological  activity  is  second  order 
with  respect  to  log  P or  n , then  log  Pg  or  nQ  can  be  calculated  and  a 
maximum  value  can  be  assigned  to  log  P or  n for  maximum  biological 
activity.  Illustrative  examples  of  this  will  be  found  in  references  18, 

21 , 23,  35,  54,  58,  61 , 66,  67,  68. 

Analgesic  and  excitant  activities  on  cats,  general  respiratory 
activities  on  rats  and  lethal  doses  on  mice  for  a number  of  morphine 
alkaloids  have  already  been  calculated  from  equations  144-147,  150,  151, 

153  and  153a.  Equations  145,  147,  151  and  153a  were  used  unsuccessfully  to 
calculate  the  respective  activities  for  six  5-alkyl dihydroworphinones , XXX 
(ref  84  page  1425).  As  the  R group  of  XXX  gets  larger,  the  deviation  between 


R=H , CH3 , C2H5,  1-C3H7, 
n-CjHn,  C6H5 


XXX 


•'mum 
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calculated  and  observed  log  A values  Increases, which  suggests  that 
steric  interaction  between  the  5-alkyloihydroroorphinone  and  the 
receptor  plays  a very  important  role  In  determining  the  biological 
activities  of  the  members  of  this  series.  Attempts  are  now  made 
in  equations  269  to  276  to  evaluate  the  steric  effects  by  E$  from  the 
data  in  Table  4.  Correlation  coefficients  for  log  Aexcit.Ca  and 
log  Ald_h  are  good  while  those  for  log  Aanalg_Ca  and  log  Agen.depr.-R 

are  not. 


lQ9  *analg-Ca  * 
^°9  Aana1g-Ca 


+ 0.25  log  Pm  + 0.41  E$  '6.05 

n = 6,  r 3 0.45  

-1.24  (log  )2  + 1.90  log  P,u  + 1 .18E$ 

-6.03 


269 


log  Aexcit-Ca 
1o9  Aexcit-Ca  2 


n * 6,  r 3 0.56  

-0.25  log  P'"  + 0.43  E$  -6.27 

n = 6,  r = 0.82  

-0.88  (log  Pm  )2  + 0.92  log  P'M  + 0.97  E$ 

-6.29 


270 

271 


1og  Agen.depr.-R 
109  Agen.depr.-R 


n 3 6,  r 3 0.85  

-0.18  log  P"  + 0.36  E$  -5.16 

n 3 6 , r 3 0.61 

-1.40  (log  P"‘  )2  + 1-59  log  P,n  + 1.89  E$ 

-5.20 

n = 6,  r 3 0.70  


Ill 

111 

274 


log  *,n  m a -0'40  l09  P“+  °'Z0  ES  -4-0’ 

LD'M  „ - 6.  r > 0.98 


275 


log  ft,  „ u - -0.17  (log  P"  V -0.18  log  P-  * 0.31  E$  -4.01 
LD'M  „ . 6,  r - 0.98 


The  sign  of  log  P"‘  terns  In  equations  271.  273  and  275 
is  negative  so  greater  lipophilicity  leads  to  greater  activity. 


The  values  for  Er  were  taken  from  reference  44. 
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The  (tog  Pm  )2  term?  in  equations  270,  272  and  274  reveal  tnat  Increase 
in  log  fJ,,‘  in  a positive  sense  contributes  to  a larger  log  A term  in  a 
negative  sense  to  a limit* and  further  increase  leads  to  > decrease  in 
activity  (parabolic  function).  The  positive  coefficient  of  Eg  indicates 
that  as  Eg  becomes  larger  in  a negative  sense  so  does  the  log  A 
term  (more  active).  Increase  in  size  of  the  R group  of  XXX  favours 
increased  activity.  Introduction  of  an  (£<.)2  term  into  equations  269p 
271  and  273  leads  to  equations  269a,  271a  and  273a. 


log  Aan  , - to. 013  log  Pn+1.?5  (E  J2  -1-22  Ec  -S.51 

3*10  "*^c.  o j 


n = 6 , r = 0.94  269a 

109  Aexcit-Ca  " “°-37  1c3  P"  +C-89  (Es>2  "0>4°  ES  *6-50 

n - 6 , r = 0. 93  271a 

'“9  V-depr-s  ’ -°  11  '°9  p"  ■°-‘17  <ES)7  *0'79  ES  -5'04 

n = 6,  r = n.65  273a 


Steric  factors  do  not  contribute  materially  to  the  degree  of  activity 
in  log  A, but  they  do  in  determining  the  magnitude  of  log  Aana*g_Ca 
ami  log  Log  A,nalg_C5  and  log  Aexcit.u  are  parabolic 

functions  of  both  Eg  and  log  P'*.  As  might  be  expected,  then,  for 
equation  269b  whore  (log  P"‘  )2  and  (Eg)2  terms  are  inserted  into 
equation  269, the  correlation  should  be  better  than  that  for  equation 
269a.  The  correlation  coefficient  of  0.45  for  equation  269  has  now 

Tog  Aana1g_Ca  = -1.13  (log  p'"  )2  + 1.51  log  P + 1.71  (Eg)2 

-0.48  Eg  -6.53 

n = 6,  r - 0.99  269b 


risen  to  0.99.  This  emphasizes  the  importance  of  steric  factors  in 
determining  the  extent  of  analgesia  by  these  compounds.  Increase  in 
the  size  of  the  group  in  the  upper  part  of  XXX  favours  increased 
activity  on  two  counts  (hydrophoDic  and  steric)  so  it  is  not  surprising 
to  find  chat  the  analgesic  activity  of  XXXI  (R'=n-C3H7)is  12,000  times 
greater  than  that  of  morphine  (193).  There  are  two  values  for  Eg  for 
the  phenyl  group.  The  value  of  +0.23  refers  to  the  thickness  of  the 
phenyl  group  while  -2.58  refers  to  the  breadth  of  the  benzene  ring. 
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The  value  of  -2.58  gives  completely  irrational  results  so  this  gives  some 
insight  into  the  space  requirements  at  the  receptor  sites. 

From  equation  272,  log  PQ“  for  this  series  of  5-alkyldihydro- 
rrorphinones  Is  <-0.52.  All  other  things  being  equal,  then,  maximum  activity 
should  be  attained  by  the  i-propyl  homolog,  but  because  of  the  large  positive 
coefficient  of  E$  in  equation  271  this  factor  exerts  a pronounced  Influence 
in  determining  the  magnitude  of  log  ^exc^t-Ca‘  Po"  ^or  ^ose 

these  alkaloids  on  mice  is  -0.53. 

Mapping  of  Enzyme  Surfaces 

The  early  model  of  drug-enzyme  or  drug-membrane  interaction,  namely 
the  "lock  and  key”  concept,  served  a useful  role;  however,  recent  developments 
have  shown  that  protein  structure  and  conformation  is  governed  to  a large 
extent  by  weak  bonds  such  as  hydrophobic  bonding  aided  by  hydrogen  bonding 
and/or  dipolar  interactions.  As  a result  the  rigid  "lock  and  key"  concept 
has  given  way  to  one  where  enzymes  have  sites  involving  greater  fluidity. 

This  is  supported  by  the  fact  that,  at  least  over  limited  ranges,  there  is 
a linear  relationship  between  Interaction  and  log  P (equations  AO,  41  anri 
211)  and  Es  (equation  46; . Such  a relationship  would  hardly  be  expecteo  on 
the  "lock  and  key”  hypothesis.  The  42  compounds  embodied  in  the  data  for 
equation  211  include  simple  alcohols,  phenols  and  anilines  as  well  as  bulky 
molecules  such  as  camphorquinone , XXXII,  neopentyl  alcohol  and  hydroxy- 
adamantane,  XXXIII. 


XXXII 


XXXIII 
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The  fact  that  rigid  bulky  molecules  fit  .nation  211  as  well  as  flexible 
compounds  argues  against  binding  by  l '<-*  and  key  interaction.  While 
some  conclusions  on  enzyme  space  were  '’r  .wn  prior  to  1S7Q,  the  majority 
of  the  work  has  appeared  since  1973. 

period  1966  - 1968 

Equations  277  and  278  were  developed  by  Hansch  at  al  (37,15) 
from  the  data  of  Metcalf  and  Fukuto  (125)  for  the  substituent  effects 
of  R on  the  inhibition  of  cholinesterase  by  alkylphosphonic acid  esters, 
XXXIV.  Comparison  of  equation  278  with  277  reveals  that  neither 


OC?  H c 


Log  K - 3.74  E$  + 7.54 

n * 13,  r = 0.901  277 

log  K - 0.15  n -1.68  a*  + 4.05  E$  + 7.21 

n » 13,  r = 0 307  278 


hydrophobic  (It)  nor  electronic  (a*)  factors  operate  in  the  inhibitory 
process.  Thus  it  would  appear  that  the  phosphonate  esters  interact 
with  the  enzyme  in  such  a way  that  contact  is  not  made  with  hydro- 
phobic regions  of  the  enzyme. 

Early  in  their  classical  work  on  the  inhibition  of 
dl hydrofolate  reductase  by  pyrimidines,  Baker  and  Shapiro  (194) 
concluded  that  inhibitors  are  bound  to  the  enzyme  by  interaction  of 
ring  electrons  with  an  electron-deficient  site  and  by  hydrophobic 
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interaction  of  one  or  more  sloe  chains.  Equations  279  to  iai  support  these 


conclusions  In  a quantitative  way  (17). 
log  jr  - -5.162  0 -5.002 

n 3 16,  r * 0.760  279 

log  J-  » 0.302  n -1 .970 

n = 16,  r = 0.328  280 

log  J-  « 0.457  r,  -5,820  a -6.951 

n = 16,  r * 0.903  281 


Introduction  of  n*  a2  or  n2  and  a2  terms  into  these  equations  made  no 
significant  improvement  in  the  correlation.  The  positive  coefficient  of  n 
in  equations  280  and  281  reveals  that  the  more  lipophilic  the  side  chain 
the  more  effective  the  inhibitor.  In  like  vein,  the  negative  coefficient 
of  a reauires  that  the  more  the  substituent  releases  electrons  to  the  ring 
the  more  effective  the  derivative  is  as  an  inhibitor.  This  is  in  accord 
with  the  conclusions  of  Baker  and  Shapiro. 

By  density  gradient  centrifugation,  Fouts  (195,196)  separated 
the  smooth-surfaced  (s)  particles  from  the  denser  -ough-surfaced  (r) 
particles  of  the  endoplasmic  reticulum  of  liver  cells  and  determined  the 
rate  of  metabolism  of  various  drugs  on  both  particles.*  Fouts  concluded 
that  the  difference  in  rates  of  metabolism  by  the  two  particles  was  due 
to  different  concentrations  of  the  enzymes  in  the  two  particles. 

Lien  and  Hansen  (25)  developed  equations  282  and  283  for  the 
ratio  R(r/a)  of  enzyme  activity  of  the  two  types  of  particles  from  rabbit 
liver  microsomes  prepared  by  the  Rothschild  and  the  Dallner  metlods. 


* 


There  was  a difference  in  the  rate  of  metabolism  of  drugs  on  the  two 
particles  from  one  to  ten  depending  on  the  drug. 


log  R,r/$)  * *0.101  log  P + 0.859 

n - 10,  r • 0.885,  s - 0.116 

log  R(r/S)  * *0.091  log  P * 0.553 

n - 7,  r * 0.923,  s * G.089 
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282 

283 


The  above  equations  permit  an  alternative  interpretation.  The  same  concen- 
tration or  activity  of  the  enzymes  can  be  present  in  both  particles  but  the 
availability  of  the  drug  to  the  enzyme  is  rate  limiting,  and  depends  upon  the 
environment  in  which  the  enzyme  is  set.  From  equation  282,  the  more  lipophilic 
the  drug  the  smaller  the  ratio  in  the  rate  of  metabolism  by  the  two  particles. 
If  some  of  the  enzymes  in  each  particle  are  in  a very  lipophilic  surrounding 
this  would  explain  why  lipophilic  drugs  like  benzpyrene  are  metabolized  at  a 
ratio  near  one.  On  this  basis  the  smooth  particles  must  have  as  well  a 
set  of  enzymes  in  a hydrophilic  medium  which  would  account  for  the  high 
ratio  for  drugs  with  a low  log  P value. 


period  1974  to  1976 

In  1974  Hansch  and  Silipo  (58)  turned  their  attention  to  Baker’s 
expanding  data  on  diami no-1 ,3,5-triazines , XVII,  as  inhibitors  of  dihydro- 
folate reductase.  Eighty-three  derivatives  of  XVII  were  included  In  the 
analysis  with  substituents,  X*,  at  C3,  . C4  and  C3C4.  Early  in  ths  work  it 
became  apparent  that  there  was  a large  difference  In  influence  of  hydrophobic 
properties  of  substituents  at  C3  and  C4  as  was  the  case  for  bulk  evaluation 
by  MR3  and  MR4 . Electronic  effects  evaluated  by  la  contributed  nothing  to 
improving  the  correlation.  Linear  combination  of  10  parameters  [n3,  n4  , 
"3,1,,  HR3,  W4,  MR3)4,a3i4,  (u3  )2,  (n4)2  and  (n3>4  }2]produced  783 
equations.  The  equation  with  two  variables  having  the  lowest  standard 
deviation  was  284  vrtiile  the  analogous  three  parameter  equation  was  285. 


* These  substituents  ranged  in  size  from  H to 


— 0CN. 


SO-O 


\=/ 
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log  Jr  = 0,464  n3+  0.181  MR4  + 6.613 

n = 83,  r x 0.834,  s = 0.422  .... 

log  l = .0.127  (n3)2  + Q.890  n3  + 0 .150  MR4  + 6.618 

n = 83,  r = 0.905,  s = 0.328  ....  285 

Equation  285  accounts  for  82%  of  the  variance  in  log  £•.  Failure  to  account 
for  18S  is  not  surprising  considering  the  enormous  variation  in  the  size  of 
substituents,  X,  in  XVII. 

Since  a plays  no  significant  role  in  the  equations,  the  electronic 
nature  of  the  substituents,  X,  need  not  be  considered.  Furthermore  n4  and 
MR3  do  not  contribute  significantly  to  any  of  the  equations.  The  reason 
why  n4  does  not  contribute  materially  to  log  ^ in  these  equations  is  that 
n (1-octanol -water)  is  not  a good  model  for  the  hydrophobic  bonding  of  the 
C4  substituent  (see  page  16).  From  equations  284  and  285  Hansch  (58) 
suggests  that,  “there  are  two  kinds  of  substituent  space  (meta  and  para)  in 
or  on  the  enzyme.  Functions  in  the  3 position  appear  to  be  placed  in  a 
typical  hydrophobic  milieu.  The  coefficient  of  about  1 for  this  term  is 
observed  quite  commonly.  Substituents  in  the  4 position  appear  to  be  thrust 
into  a more  apolar  region  which  n4  does  not  model  well.  It  seems  likely 
that  groups  in  the  4 position  cause  inhibition  by  producing  conformational 
changes  in  di hydrofolate  reductase  by  more  firmiy  attaching  the  inhibitor 
through  dispersion  forces  or  by  a combination  of  both."  Adding  an  (MR4)2 
term  to  the  equations  does  not  improve  the  correlation  so  log  ~ is  not  a 
parabolic  function  of  MR ,so  still  bulkier  substituents  could  be  accom- 
modated where  C4~substi tuents  interact  with  the  enzyme.  On  the  other 
hand  (n3)2  terms  contribute  materially  to  equation  285  so  log  J-  is  a para- 
bolic function  of  li  and  it  is  unwise  to  introduce  a group  at  C3  with  a n 
value  greater  than  R0  whi,.h  is  +3.5.  To  examine  the  effect  of  bulkier 
groups  at  C4  in  the  di ami  no-1 ,3 ,5-tri azi nes  , Si  1 i po  and  Hansch  (69)  extended 
this  study  to  include  the  analysis  of  the  interaction  of  2d4  of  these  com- 
pounds with  dihydrofolate  reductase. 

Sill po  and  Hansch  (66)  have  analyzed  the  data  of  Baker  and 
Wood  (197,198)  on  the  inhibition,  Jr,  of  xanthine  oxidase  by  thirty 
9- (substi tuted-phenvl ) guanines,  XXXV. 


XXXV 
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These  workers  developed  511  equations  involving  8 variables  chosen 
from  n2,  \\ 3 , r.4,  n3j4>  in,  (in)2  MR2 , MR3 , MR4,  MR3j4» 

(MR3|4)2  IMR,  (iMR)2,  E$  , Es_4  , a,  o+,  £3,  f 4 , I?,  Jfia  . . 

0 and  D2 . In  all,  more  than  2,000  equations  were  examined.  The 
equation  with  the  highest  correlation  coefficient  and  lowest  standard 
deviation  is  equation  286.  The  relative  importance  of  the  various 
terms  in  equation  286  is  manifst  in  the  statistical  data  for  equations 
287  to  290. 


log  1 = 0.203  MR3>4  + 1.259  E$_2 

log  J-  ■ 0.24  E-_„  + 6.12 

log  i = 0.26  MR3>4  + 5.90 

log  i - 1.13  Es_2  + 4.98 

lo9  1.47  Es_2  * 0.40  Es_„ 


+ 0.432  Es_4+4.327 


n - 

30, 

r = 

0.924, 

s = 

0.228  . 

286 

n = 

3 0 ; 

r = 

0.359, 

s = 

0.535  . 

287 

n = 

30, 

r = 

0.476, 

s = 

0.505  . 

288 

n = 

30, 

r = 

0.630, 

s = 

0.445  . 

289 

n = 

30. 

r = 

0.857, 

5 = 

0.301  . 

290 

The  outstanding  observation  to  be  made  from  equation  286  is  that  the  9-phenyl 

ring  of  these  compounds  does  not  locate  itself  in  hydrophobic  space  on  the 

enzyme.  In  no  case  did  n give  as  good  results  as  MR  (see  col  linearity 

page  95).  Equation  286  accommodates  large  substituents  (e.g.f  3-NHC0C6H5  and 

4-0CH2CH2CH2C6H5)  in  the  meta  and  para  position  indicating  the  great  flexibility 

in  the  enzyme  space  around  these  positions.  Since  inclusion  of  an  (MR3  u)z 

\ ’ 

term  in  the  equation  did  not  improve  the  correlation,  then  log  ^ does  not 
appear  to  be  a parabolic  function  of  MR3  4 so  still  more  active  derivatives 
can  be  prepared  by  using  still  larger  substituents  in  the  meta  and  para 
positions.  The  coefficients  of  E<-_2  and  Ej._4  in  equation  286  are  positive 
and  this,  combined  with  the  fact  that  the  larger  the  group  the  larger  E<- 
values  are  in  a negative  sense,  implies  that  there  are  obstacles  in  the 
C2  and  C4  space  on  the  enzyme  which  hinder  the  entry  of  large  groups  in 
these  positions  in  XXXV.  The  role  of  MR  for  substituents  is  an  ambivalent 
one.  MR  may  be  evaluating  dispersion  forces  binding  the  inhibitor  to  the 
enzyme  (198)  or  it  may  gauge  the  volume  of  the  substituent  and  its  ability 
to  distort  the  conformation  of  the  enzyme  so  as  to  preclude  interaction 
with  the  proper  groups  on  the  inhibitor. 

Inhibition  of  guinea  pig  complement*by  benzyl pyri di ni un  ions 
(199-201)  and  benzamedines  (202-206)  has  been  studied  by  Baker  and  his 
collaborators  and  these  results  have  been  analyzed  by  Hansch  et  al  (60,67). 


* Complement,  as  described  by  Hansch  and  Yoshimoto  (60),  "consists  of 
11  distinct  proteins  which  are  required  for  cell  lysis  brought  about 
via  antibodies  and  complement.  The  function  of  the  antibody  is  to 
identify  the  invading  cell  as  a foreign  organism  and  activate  complement 
attack  which  results  in  cell  lysis  by  means  of  the  proteolytic  enzymes." 
This  mechanism  operates  in  the  rejection  of  foreign  tissue  and  organ 
transplants. 


few 
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From  the  complex  nature  of  the  benzamidines,  XXXVI,  it 
was  apparent  that  durnny  variables  would  have  to  be  employed  to  encompass 
the  diverse  structures  present  in  the  108  benzamidines. 


XXXVI 


The  structural  features  of  substituent  X in  XXXVI  were  para  meterized 
by  Hi  , C],  and  MRj  while  n2,  o2  and  MR2  were  assigned  to  the 
Y substituent  in  the  other  benzene  ring.  The  indicator  variable 
Di  characterized  the  grouping  inside  the  box  of  XXXVI.  For  XXXVI 

compounds  with  ■?  bridge  L-0(CH2 )20- , -(CH2)4  -0(CH2)30-,  -0 ( CH2 ) 3 - » 

-0(CH2  )40- , -0(CH2 )it-l  to  a second  ring  I I*D3  was  set  at  1.00.  For  all 

other  inhibitors  was  assigned  a value  of  0.00.  D2  was  introduced  to 

categorize  a pyridine  ring  at  the  end  of  the  side  chain.  The  connecting 

side  chain  was  evaluated  by  n or  MR.  There  was  an  abnormal  exaltation 

in  inhibitory  activity  when  the  grouping  | was  attached 

NHC  - Z - C&H5 

at  C3  of  ring  1 t nf  XXXVI.  When  7.  was  zero,  NH,  CH?  , NHCH2CH?  or 

CH20,  D3  was  assigned  a value  of  1.00.  When  these  groups  are  at  C2  or  C4 

of  ring  Htheir  contribution  was  accounted  for  by  n and  MR. 


* includes  the  terminal  benzene  ring, II. 


108 


This  analysis  (60)  was  approached  in  the  same  way  as  outlined  above, 
and  briefly  it  may  be  stated  that  equations  291  and  292  provided  the  best 
fit  for  the  data.  The  correlations  from  equations  291  and  292  are  so  similar 

log  ~ = 0.146  MRj  f2  + 1 -068  Dj  + 0.520  02  * 0.429  05  + 2.425 

n --  108,  = 0.935,  s = 0.258.  . . 291 

log  - = 0.211  nli2  + 1 .345  Dj  + 0.620  D2  + 0.565  D3  + 2.440 

n = 108,  r = 0.931  , s = 0.267.  . . 292 

that  it  is  not  possible  to  say  with  certainty  whether  the  substituent 
effect*  is  hydrophobic  in  nature  or  due  to  the  polarizability  of  the  sub- 
stituents. Hansch  and  Yoshimoto  (60)  as  well  as  Coats  (207)  present  circum- 
stantial evidence  favouring  MRlt2  as  the  descriptor  and  not  n.  Expansion 
of  equations  291  and  292  to  include  (KRlt2  )2  or  (nll2)2  terms  did  not  reduce 

the  variance  in  log  g-  so  log  ^ is  not  a parabolic  function  of  MR1>2  or  n1>2- 

The  positive  coefficient  of  HP.j  2 indicates  that  the  larger  the  substituents 
X and  Y,  the  more  effective  is  the  inhibitor.  Resolution  of  MRi>2  into  MR} 
and  MR 2 gave  an  equation  in  which  the  coefficients  of  these  two  terms  were 
the  same  sign  and  about  the  same  magnitude,  so  the  inhibitory  effect  due  to 
the  substituents  in  each  ring  is  about  the  same.  Hansch  and  Yoshimoto  con- 
cluded that  MR  ;s  not  reflecting  conformational  changes  in  complement  by 
the  substituents  but  is,  more  likely,  evaluating  the  binding  of  the 
inhibitor  to  the  complement  by  dispersion  forces.  Hansch  and  Yoshimoto 
then  proceeded  to  analyze  the  dummy  parameters  Dj , D2  and  D3. 

The  inhibition  of  complement  by  69  benzyl pyridini urn  ions  (67) 
has  been  treated  similarly  as  has  that  (74)  of  chymotrypsin , trypsin, 
thymidine,  phosphoryl ase , uridine  phosphoryl ase, thymidilate  synthetase, 
cytosine,  nucleoside  deaminase,  malate  dehydrogenase,  glutamate  dehydro- 
genase, lactate  dehydrogenase  and  glyceral dehyde-phosphate  dehydrogenase. 

The  QSAR  for  the  inhibition  of  malate  and  glutamate  dehydrogenase  by 
1 ,4-dihydro-4-quinolone-3-carboxylate  ions,  XXXVII,  are  very  similar  as 
seen  from  equations  293  and  294. 


★ 


The  degree  of  coll  inearity  between  P.  and  MR  is  presented  in  the 
correlation  matrix  in  Table  III  of  reference  60. 


H 

XXXVII 


log  • 0.699  n5  + 0.290  MR6>7|8  -1.121  Ix  + 3.156 


n = 75,  r = 0.943,  s = 0.385  293 

log  jl  = 0.491  n5  + 0-233  MR6  -0.553  Ix  + 3.355 

n = 87,  r = 0.948,  s = 0.253  294 


From  the  equations  derived  for  the  inhibition  of  these  two  enzymes  by 
1 ,4-dihydroquinolone-3-carboxyl ates  Yoshimoto  and  Hansch  (74) 
presented  the  schematic  drawings  XXXVIII  and  XXXIX  for  enzymic  space 
relative  to  the  1 ,4-dihydroquinolone-3-carboxylates . A similar 


XXXVIII 


XXXIX 


no 


diagram  was  presented  (74)  for  the  much  more  complex  process  of  inhibition 
of  di hydrofolate  reductase  by  diamino-1 ,3, 5-triazines. 

Systematic  Drug  Design 

The  iaeal  chemotherapeutic  agent  is  one  that  is  non-toxic  to 
humans  but  quite  toxic  to  parasitic  or  malignant  cells.  A most  rational 
approach  to  this  study  is  the  examination  in  Isolation  of  enzymic  or 
membrane-controlled  processes  in  both  host  and  parasite.  If  a chemo- 
therapeutic is  known  for  these  processes  then  systematic  drug  modifi- 
cation can  start  with  this  compound  or  with  one  with  similar  chemical 
properties  and  log  P value.  The  systematic  modification  will  involve 
planned  molecular  changes  that  become  progressively  more  extensive.  The 
extensive  and  classical  work  of  Baker  and  his  collaborators  illustrates 
this  approach  which  he  has  subdivided  into  four  steps  which  have  been 
concisely  suumarized  by  Hansch  (40). 

1)  "An  enzyme  must  be  selected  and  a reversible  inhibitor*  found. 
Modification  of  its  chemically  more  active  groups  will  lead  to 
suitable  reversible  inhibitors.  Binding  points  on  the  reversible 
inhibitor  that  complex  with  the  enzyme  should  be  determined. 
Careful  systematic  variation  of  substituents  yields  the  necessary 
insight. 

2)  Areas  on  the  inhibitor  should  be  determined  in  which  bulky  groups 
can  be  placed.  This  uncovers  two  types  of  positions: 

a)  large  flexible  hydrophobic  areas  and  b)  non-contact  areas 
between  inhibitor  and  enzyme. 

3)  Once  the  non-contact  area  is  determined,  then  a group  that  can 
form  a covalent  bond  with  common  enzymic  functions  should  be 
placed  in  this  area.  The  length  of  the  side  chain  by  which  the 
function  is  attached  to  the  parent  inhibitor  can  be  varied  so 
that  the  function  can  react  irreversibly  with  a group  on  the 
enzyme  outside  the  active  site. 


4)  After  finding  the  ideal  length  and  flexibility  of  the  side  chain 
which  is  to  react  irreversibly,  then  variations  in  the  active 
function  itself  should  be  investigated  in  order  to  find  the 
function  with  the  ideal  stereoelectronic  specificity,  i.e., 
one  whose  reactions  with  the  other  molecules  in  the  host  are 
minimized, " 

Careful  systematic  variation  of  substituents  has  been 
approached  in  two  ways.  1)  Topliss  (96,97)  has  developed  a systematic 
approach  to  the  selection  of  six  substituents  that  will  reveal  the  sub- 
stituent trend  as  related  to  biological  activity.  2)  Single  parameter 
equations  developed  from  few  well  selected  derivatives  will  reveal  the 
influence  of  substituents  upon  biological  activity.  The  application 
of  the  Topliss  tree  to  the  study  of  substituent  trends  is  clearly  demon- 
strated by  Hansch's  application  to  the  study  of  ellipticene  derivatives 
based  upon  the  parent  compound,  XL.  (41). 


Move  1 in  Fig.  IX  is  to  introduce  a bromine  atom  at  C5  in  XL  and  deter- 
mine the  biological  activity  of  the  5-bromo-deri vati ve . The  5-bromo- 
derivative  can  be  more  active  (+),  less  active  (-)  or  about  the  same 
activity  (0)  as  the  parent  compound.  If  the  5-bromo-deri vati ve  is  more 
active  than  XL  (right  hand  side  of  Fig.  IX),  then  the  5-SCF 3*  deriva- 
tive of  XL  is  prepared  and  tested.  If  it  is  more  active  than  the  5-SCF 3 
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deri vati ve»movement  is  to  the  left  and  the  5-phenyl  derivative  is 
prepared;  if  less  active .then  the  5-CH3  derivative  is  prepared. 
Introduction  of  the  phenyl  group  enhances  r:  but  not  a relative  to 
5-SCF3.  The  n and  a constants  for  5-CH3  are  lower,  in  a positive 
sense,  than  are  those  for  5-SCF3. 

If  the  biological  activity  cf  the  5-bromo-deri vati ve  is  less 
(-)  than  that  of  XL,  movement  in  Fig.  IX  is  to  the  left  and  the  5- 
hydroxy-deri vati ve  of  XL  is  prepared  and  tested.  If  the  5-hydroxy- 
derivative  is  more  active  (+)  then  a more  water  soluble  one  still,  the 
5-NH2CQNH-de,'i  vati  ve  could  be  prepared  and  tested,  if  the  result  of  the 
5-OH  derivative  is  (-),  it  could  De  that  a compound  with  a positive  0 
and  a negative  n value  is  needed.  The  5-CN  would  meet  these  requirements 
nicely, 

In  the  study  of  the  inhibition  of  complement  by  benzami dines , 
Hansch  and  Yoshimoto  (53)  recommend  the  early  attempt  to  establish 
thermodynamic  relationships  by  development  of  single  parameter  equations 
on  a few  well  chosen  compounds.  Equations  of  the  form  of  equations 
295-297  could  have  been  developed  long  before  108  compounds  had  been 


synthesized. 

log  J-  = 1 .28  Di  * 2.80 

n = 108,  r = 0.785,  s « 0.447  295 

log  ~ = 1 .05  Dj  + 0.19  MRL  + 2.42 

n = 108,  r = 0.905,  s = 0.308  296 

log  = 0.99  D,  + 0.16  MR1>2  0.41  D3  + 2.48 

n = 108,  r = 0.927  , s = 0.273  297 


This  would  have  led  to  better  derivatives  and  eliminated  redundancy  in 
the  synthesis  of  derivatives.  As  a result, equations  291  and  292  would 
have  been  obtained  in  less  time  and  with  less  work. 

Methods  for  mapping  enzymes  relative  to  the  inhibitor  have  been 
discussed  on  pages  100  to  110. 

The  concurrent  examination  of  enzyme  or  membrane  in  both 
host  and  parasite  will  ensure  attainment  cf  the  highest  degree  of 
selective  toxicity.  Realization  of  a highly  active  inhibitor  on 
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isolated  enzyme  does  not  ensure  high  activity  in  in  vivo  tests  (see 
page  84).  Competing  processes  such  as  metabolism,  elimination  or 
anchoring  of  the  inhibitor  on  the  way  to  the  active  site  may  reduce  the 
efficacy  of  the  drug  in  in  vivo  experiments  tc  almost  zero.  Log  Pq  for 
another  series  of  drugs  acting  at  the  same  site  in  the  same  animal  or 
that  of  the  same  series  of  drugs  acting  analogously  in  another  animal 
(see  pages  86  to  90)  provides  a reasonably  safe  guideline  to  follow  in 
selecting  the  member  of  the  series  for  the  most  favourable  “random  walku 
in  the  in  vivo  experiments. 

In  general,  high  rates  of  metabolism*  and  elimination  are  favoured 
respectively  by  highly  lipophilic  and  hydrophobic  compounds.  Equations 
can  be  developed  for  these  processes  and  the  sign  and  magnitude  of  the 
coefficients  of  log  P ccr  3ared  with  log  PQ  for  this  series  of  compounds 
on  this  in  vivo  biological  process  (see  pages  83  to  90).  Employing  these 
methods  led  Hansch  and  Silipo  (58)  to  the  synthesis  of  the  diamino-1 ,3,5- 
triazine,  XLI,  as  a more  effective  inhibitor  of  dihydrofolate  reductase. 


* Metabolites,  in  general,  are  quite  inactive  but  compounds  like  the  tetra 
hydrocannabinol s (192)  can  be  encountered  where  the  primary  metabolite  i 
the  active  constituent. 
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In  conclusion,  methods  have  been  developed  for  the  calculation 
of  the  biological  activities  of  a series  of  drugs  cn  one  organism  from 
their  observed  activities  on  another  organism  (see  pages  69  to  73  ). 
Preliminary  work  (84)  indicates  that  this  method  can  be  applied  to  the 
calculation  of  biological  activities  of  a series  of  drugs  on  man  from 
their  observed  activities  on  test  animals. 
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Chymo trypsin 
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- see  collinearity 
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reductase  by,  103 
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Dihydrofolate  reductase 
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- morphine,  55 
Dipole  moments , 61 
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- mapping  of,  100 


Equations  for 
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- 109  Aanalg-H*  67 
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reflex,  56 
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- effect  on  log  P,  18 
Hydrophobic  bonds,  18,  21,  82 
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- specific  and  non-specific,  81 
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dehydrogenase,  108 
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- equations  for  inhibition  of 
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- equations  for  inhibition  of 

growth,  89 

Guinea  pig  complement 
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k , 56 
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krf,  10,  53.  60 
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Lactate  dehydrogenase 

- inhibition  of,  108 
Lock  and  key,  100 
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- equations  for  5-al kyldi hydro- 

morphinones,  97 

- equations  for  morphines,  57 

103  Aanalg-M 

- equations  for  morphines,  67 

Log  abm 

- equations  for.  74 
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- equations  for  5-al kyldi hydro- 

morphinones,  97 

- equations  for  morphines,  68 
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- equations  for  5-al kyldi hydro- 

morphinones,  97 

- equations  for  morphines,  68 

Log  \onv-M 

- equation  for  morphines,  68 
Log  ALD-M 

- equations  for  5-al kyldi hydro- 
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- equations  for  morphines,  58 
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- equations  for  morphines,  68 
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- equations  for,  56,  58 
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- equations  for,  60 
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- equations  for  S.  aureus,  62 

Log  IGso+a 
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- dimeric  oxidation  of,  55 
Morphine  alkaloids 

- biological  test  methods  for,  65 
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for  analysis,  67 
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Non-barbiturates 
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- monamine  oxidase  inhibitors,  47 
9-Phenylguanlnes 

- inhibition  of  xanthine 
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- antlroalarial  activity  of,  80 
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- enzymic  hydrolysis  of,  37 
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- inhibition  of,  108 
I'.,  8,  12,  15.  18,  105 
nQ,  83,  % 
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Paeudomonaa  aeruginoaa , 64 
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Pyrimidines 

- inhibition  of  dihydrofolate 
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- appl icatlons  of,  95 
QSAR  equations 

- from  de  novo  substituent 

constants,  80 
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considerations,  23 

- from  experimentally  determined 

physico-chemical  constants, 
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Rabbit  kidney  reductase,  37 
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specific,  65 
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Cj , 15 

aR  * 15 
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Symbols,  3 
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- inhibition  of,  108 
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Topi i ss  tree , 111 
Toxi city 
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